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CHAPTEI;  I 

' • SELF  DIPTUSION  Hi  IMPURE  SODIUM  CHLORIDE 

Ao  Introduction 

Because  of  the  simplicity  of  their  crystallographic 
structure , the  alkali  halides  occupy  a prominent  position  in 
I the  field  of  solid  state  physicso  X-ray  diffraction  experiments 

and  all  supplementary  evidence  Indicate  that  an  alkali  halide 
j ^ such  as  sodium  chloride  has  a cubic  structure  with  sodium  and 

j f chlorine  ions  so  arranged  that  each  soditun  has  six  nearest 

I • 1 

i neighbor  chlorine  ions  and  vice  versa©  A two  dimensional 

representation  of  such  a lattice  is  shown  in  Figure  lA. 

, The  electrical  conductivity  of  the  alkali  halides  has 

. ' been  investigated  by  many  workers,  the  pioneer  work  having  been 

done  by  Lehfeldt^^^  and  Tubandt©^^^  The  latter  also  established 
the  fact  that  the  conductivity  in  alkali  halides  is  ionic  in 
I nature©  The  observed  conductivity  can  be  explained  by  assuming 

the  presence  of  imperfections  in  the  lattice.  This  was  shown 
by  the  theoretical  investigations  of  Frenkel, Wagner, 

Schottky  ,^5)  Jost,^^*^  and  later  Mott  and  Littleton©  Frenkel^^^ 
proposed  a defect  in  which  a number  of  ions  have  been  removed 
from  their  normal  sites  and  displaced  to  interstitial  positions 
in  the  lattice©  This  is  shown  in  Figure  IB,  Conduction  can 
then  take  place  by  the  motion  of  vacancies  created  when  an  ion  is 
removed  to  an  interstitial  position,  by  the  motion  of  interstitial 
. i ■ ions,  or  by  both©  A different  type  of  defect,  Illustrated  in 

I Figure  10,  has  been  proposed  by  Schottky©^^^  It  consists  of 


( 

) 


vacancies  created  in  the  lattice  by  the  removal  of  an  equal 
ntimber  of  positive  and  negative  ions  from  the  interior  to  the 
surface  of  the  crystal*  Conduction  in  this  case  can  take  place 

either  by  the  motion  of  positive  or  negative  ion  vacancies  or 

(c ) ) (7 ) 

both.  The  work  of  Schottky,'^'  Jos'-,'®'  and  Mott  and  Littleton 

has  shown  that  the  predominant  defect  in  the  alkali  halides  is 

of  the  Schottky  type* 

Beth  Frenkel  and  Schottky  defects  can  be  caused  by 
thermal  agitation,  by  the  presence  of  impurity  ions  in  the  lattice, 
or*  by  departure  from  exact  stoichiometric  proportions.  Theoretical 
expressions  for  the  concentration  of  these  defects  can  be  derived 
by  minimizing  the  free  energy. 

It  was  further  observed  by  Lehfeldt,^^^  and  later  by 

/q ) 

’loch  and  Wagner,'^'  that  when  small  amounts  of  a divalent  salt 
were  addea  to  a monovalent  alkali  halide  lattice  the  conductivity 
increased  considerably.  The  mechanism  causing  this  increase  can 
be  explained  approximately  in  the  following  way.  If  a solid 
solution  is  formed  in  which  a divalent  positive  ion  replaces  a 
monovalent  positive  ion  in  the  lattice,  then  a second  monovalent 
ion  must  leave  the  lattice  to  preserve  electrical  neutrality. 

Thus  a positive  ion  vacancy  is  created  for  each  divalent  ion 
introduced  substltutionally  into  the  lattice.  This  situation  is 
illustrated  in  Figure  ID.  In  recent  years  extensive  studies  of 
conductivity  as  a function  of  impurity  contesnt  have  been  made 
by  Etzel  and  Maurer and  by  Bean  and  Maurer. 


Another  process  occurring  in  alkali  halides  vdiich  is 


3 


due  to  the  presence  of  vacancies  is  that  of  diffusion.  The 
ions  in  the  crystal  can  diffuse  from  one  part  of  the  crystal  to 
another  by  jumping  into  a vacancy,  in  turn  leaving  a vacancy 

behind,  Self-dif fusion  of  the  cation  has  been  measured  by 

(12)  (13) 

Mapother-.  Crooks,  and  Maurer  for  NaCl  and  NaBr  and  by  V/itt 

for  KClo 

In  view  of  the  fact  that  the  same  mechanism  is 
responsible  for  both  conduction  and  diffusion,  a relationship 
between  the  two  exists.  This  was  first  derived  by  Einstein 
(see  also  appendix  I),  and  can  be  stated  as 


r 

"IT 


ne" 


(1) 


vrhere 

^ = conductivity 
D » diffusion  coefficient 
n » density  of  carriers 
e = charge  of  carriers 
k “ Boltzmann  constant 
T = absolute  temperature. 

This  equation,  however,  is  valid  only  if  both  conductivity 
and  diffusion  are  brou^t  about  by  exactly  the  same  process,  be 
it  motion  of  the  ions  or  motion  of  the  vacancies.  It  should  be 
pointed  out  that  there  are  processes  whereby  diffusion  can  take 
place  without  the  transport  of  electric  charge,  thus  invalidating 
the  Einstein  relation.  Such  processes  include  the  formation  of 
positive  ion  vacancy-negative  ion  vacancy  pairs  (Figure  IE)  and 
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divalent  ion-vacancy  complexes  (Figure  IF).  The  reason  such 
pairs  and  complexes  Invalidate  the  Einstein  relation  is  due  to 
the  fact  that  they  form  neutral  entities  in  the  lattice.  Thus 
the  positive  ion  vacancy  in  each  case  will  not  contribute  to 
the  conductivity,  but  may  contribute  to  the  diffusion.  (The 
contribution  of  the  negative  ion  vacancy  can  usually  be  neglected^ 
since  its  mobility  is  small  compar»>rt  t.r,  of  the  positive  ion 

vacancy)® 

The  validity  of  the  Einstein  relation  has  been  tested 
for  NaCl  and  NaBr  crystals  by  the  work  of  Mapother,  Crooks,  and 
Maurer. These  investigators  found  the  relation  to  be  valid 
in  the  high  temperature  or  intrinsic  region  (above  approximately 
500®C),  but  noticed  a distinct  deviation  in  the  low  temperature 
or  structure  sensitive  region.  The  latter  behavior  they  attributed 
to  the  existence  of  complexes,  involving  impurity  ions  present 
even  in  so-called  ”pure”  crystals.  In  the  intrinsic  region, 
on  the  other  hand,  the  thermally  generated  vacancies  outnumber 
those  introduced  by  impurities.  Hence  in  this  region  the  complexes 
become  rather  unimportant  compared  to  the  total  number  of 
vacancies,  and  the  Einstein  relation  is  satisfied. 

The  work  of  Etzel  and  I-laurer^^^^  previously  referred 
to  also  is  concerned  with  these  complexes.  It  was  found  that 
the  addition  of  Cd'*"^  ions  to  the  NaCl  lattice  increased  the 
conductivity,  but  not  quite  proportionately  to  the  amount  of  Cd 
added.  This  was  attributed  tc  a larger  number  of  complexes  being 
formed  per  Cd'^'*’  ions  added  as  the  concentration  is  increased. 

From  their  results,  Etzel  <and  Maurer  could  obtain  values  for  the 
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number  of  complexes  formed  per  Cd'*"*’  ions  added  (or,  as  it  is 
commonly  referred  to,  the  degree  of  association)  for  various 
concentrations  and  temperatureso 

The  purpose  of  the  present  experiment  is  to  obtain  a » 

further  check  on  the  degree  of  association  in  the  system  NaCl  ! 

: NaCl  + CdCl2o  This  is  achieved  by  measuring  not  only  the 

. . i 

condxictivities  of  this  system  for  various  concentrations  of  Cd  I 

ions  (and  incidentally  also  over  a wider  temperature  range  than 
previous  data),  but  also  the  self-diffusion  coefficients  of  the 
I I Na”*"  ion,  using  a radioactive  tracer  technique » The  conductivity 

and  the  diffbsion  coefficients  can  then  be  compared  by  means  of 

the  Einstein  relation„  Any  deviation  should  therefore  yield  ’ 

airect  information  about  the  degree  of  association  at  that 

} ■ 

'■  particular  temperature  and  concentration.  These  results  can 

then  be  compared  with  the  theor^^o 

^ ! 

' Bg  Theory 

i 

In  order  to  account  quantitatively  for  the  failure  of 
the  Einstein  relation  in  some  alkali  halides  with  added  divalent 
impurities,  it  is  essential  to  know  the  number  of  complexes  ^ 

^ present  as  a function  of  temperature  and  impurity  content,  ioe» 

the  degree  of  association . In  addition,  the  binding  energy  of  * 

j,  a complex  must  be  obtained„  V/e  shall  briefly  derive  the  degree 

I of  association,  using  first  a very  simple  approach  and  then  a 

j somewhat  more  detailed  one.  Then  we  shall  show  how  the  binding 

energy  can  be  arrived  at,, 

j The  basis  of  the  calculations  of  the  degree  of  ! 

! 

f 

! 


i 


6 


association  is  the  mass  action  law.  This  approach  was  first 
used  by  Stasiw  and  Teltow, ^^5)  later  by  Teltowy^^^^  Etzel  and 
Maurer, and  recently  by  Lidlard.'l?!  The  last  t«o  papers 
are  of  special  interest  in  connection  with  the  preserjt  work 
because  they  deal  with  the  same  system  used  here,  i..e. 

NaCl  + CdClgo  We  will  follow  in  the  main  the  paper  by  Lidiard, 
Sj.nce  it  contains  not  only  the  relatively  simple  approach  of 
the  previous  papers,  but  also  considers  the  problem  without 
some  of  the  simplifying  assumptions  made  by  the  earlier  workers. 
Let  us  take  an  alkali  halide  lattice  containing 
solvent  cations  and  divalent  impurity  ions.  At  any  given 
temperature  let  there  be  n complexes  present,  and  hence  N^^n 
unassociated  vacancies  and  N^^-n  unassociated  impurity  ions. 

Then  from  the  equilibrium  condition 


unassociated  impurity  ion  unassociated 


vacancy 


"T'  complex 


(2) 


the  mass  action  law  gives 


N^n 


(N^  - n)^ 


K(T), 


(3) 


It  is  to  be  noted  that  K(T)  is  independent  of  vdiich 
has  the  sigiificance  of  neglecting  the  mutual  interactions  between 
the  unassociated  vacancies  and  impurity  ions.  The  effect  of 
taking  these  interactions  into  account  will  be  treated  later. 

An  explicit  form  of  K(T)  is  obtained  by  minimiziiig 
the  free  energy  with  respect  to  n.  Again  neglecting  the  inter- 
actions of  the  unassociated  charges,  the  free  energy  is  given  by 


n 


. 1 
I t 


, I 


Fq  « -kT  In 


z”  TF  (V^2M...2sJ 

s»0 


n>! 


)j 

“ N^f  (nj-n)fi  (N^-n)/ 


ne 

er. 


(4) 


The  first  term  is  the  entropy  term,  z being  the  number 
of  distinguishable  orientations  of  a complex,  (For  the  cast  of 
NaCl,  the  vacancy  can  be  in  12  different  nearest  neighbor  positions 
in  the  cation  sublatticeo  Hence  a 12,)  The  second  term 
represents  the  energy  of  associations,  e being  the  static  dielectric 
constant  of  pure  NaGl.  and  r^^  the  nearest  neighbor  distance.  (For 
the  case  of  NaCl,  r^  =2  a,  \rtiere  a is  the  anion-cation  separation, 
and  e = 5 ,62  ) , 

Kiiiimizing  with  respect  to  n^  we  get 


N^n 


z exp 


( 


2 


‘(N^-n.)^ 

The  degree  of  associaticn  is  defined  by 


er^k® 


(5) 


n 

p “ -q 


(6) 


and  can  be  obtained  from  (5),  VJe  thus  find,  using  z “ 12, 


where 


p = 1 _ 1.  ( _1  + » '2y } 

y 


T 

y = exp  -,o  + In (240) 
2 

■p  „ ® — _ 

o erc£ 


and 


c 


N, 


(7) 

(S) 

(9) 

(10) 


I 

J 
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Tq  as  defined  by  (9)  i^  the  binding  energy  of  a complex  in  units 
of  Boltzmann's  constant o 

V/e  next  take  into  account  the  Coulomb  interactions  of 
the  unassociated  charges*  This  can  be  done  by  treating  these 
charges  as  being  equivalent  to  a solution  of  a binary  electrolyte 
in  a medium  of  dielectric  constant  e.  Such  systems  can  then  be 
handled  by  the  Debye-Hueckel  theory  of  electrolytes.  This  approach 
is  legitimate  even  though  the  charges  in  question  are  in  fixed 
lattice  positions,  rather  than  freely  moving  as  in  a true  solution, 
since  the  average  distance  between  these  charges  is  much  larger 
than  the  lattice  spacing.  Care  nust  be  taken,  however,  not  to  let 
unassociated  vacancies  and  impurity  ions  approach  closer  than  the 
next  nearest  neighbor  distancs.  If  they  approached  any  closer 
they  v/ould  then  of  course  be  classified  as  complexes. 

The  main  effect  of  applying  the  Deby-Hueckel  theory  to 
the  problem  at  hand  is  to  introduce  an  additional  t erm  to  the  free 
energy,  which  must  be  taken  into  consideration  when  minimizing  with 
respect  to  n.  The  details  of  this  term  as  well  as  the  new  form 
that  eqo  (5)  takes  are  contained  in  Lidiard's  paper  and  will  not 
be  repeated  here.  We  will  merely  quote  the  result  for  the  NaCi 
type  lattice  with  impurities,  from  which  the  degree  of  association 
can  be  calculated.  It  is 

r 


- exp 


\o  + ln(l2c) 


t1/2 


(Tp/T  (d-p  }c  j ' 


(11) 


\i  1»P  )cTo/r]^'/2  } 

This  equation  can  be  solved  for  p by  successive  approximations. 
(See  lidiai'd's  paper  for  details..) 


•>/  'i 


Before  we  can  apply  equations  (7)  and  (11/  to  compute  the 
degree  of  association  for  specific  impurity  concentrations  and 
temperatures,  however,  we  still  must  determine  T^,  i.e,  the  binding 
energy  of  a complex.  There  are  at  least  three  ways  of  doing  this. 
Probably  th«j  simplest  way  is  to  take  the  Coulomb  potential  energy 
of  tiiro  charges  of  opposite  sign  separated  by  a distance  , in 


a medium  of  dielectric  constant  e.  This  gives 

e2 

V a „ S 

e<2*  a 


(12) 


from  which  we  find  T^  =»  7000°K,  or  a binding  energy  of  .60  eV. 

In  arriving  at  this  value,  however,  we  have  completely  neglected 
all  interactions  except  the  Coulomb  interaction  at  the  nearest 
neiglibcr  distance.  It  will  be  shown  below  that  these  other  inter- 
actions are  by  no  means  negligible^. 

A second  way  of  finding  T^  is  due  to  Etzel  and  Maurer. 

By  v.Titing  the  conductivity  as 


(T  - 


. 2e^ 


-UAT  > 
e ' c ( 1-p } 


(13) 


where  V is  the  jump  frequency  of  a vacancy  arid  U the  activation 
energy  for  such  a jump,  we  find,  substituting  for  p from  (7),  that 

_ ..  To/r  , _.,2 


a(r  + 12e 


(a(T)' 


where 


(14) 


(15) 


Etzel  and  Maurer  then  fitted  their  data  of  (T  vs„  c to  a parabola 


c = F<r  + L(T" 


(16) 


1 
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By  comparing  (I4)  with  (16),  it  is  found  that; 
12  =-  L/F^ 


(17) 


from  which  can  be  found.  turns  out  to  be  3760  K,  corres- 
ponding to  a binding  energy  for  the  complex  of  .33  eV, 

There  are,  however,  some  shortcomings  in  this  method. 

First  of  all,  we  have  neglected  the  effect  of  the  atmosphere  of 
unassociated  impurity  ions  on  the  motion  of  the  vacancies  under  the 
influence  of  a field.  This  effect  will  be  a hindrance  of  the 
motion,  thus  CT  must  be  multipl.*-^d  by  a factor  (1-f),  where  f<l. 


An  expression  for  f has  been  given  by  Pitts, 


(16) 


When  applied  to 


an  NaCl  type  lattice,  it  results  in  the  following  equation  for  o'  : 

0^2  . -UAT  , . f 

^ ~ alcT"  ® c(l-p)  n 1 - 

(16) 

- V >A_  / -^  ; 


Secondly,  in  the  previous  method,  the  Interaction  between  the 
unassociated  charges  have  been  neglected.  This  can  be  corrected 
by  substituting  into  (16)  the  expression  for  p derived  from  (11). 
However,  we  then  no  longer  find  a simple  relationship  of  the  form 
(14)0  Specifically,  a plot  of  vs.  ^ gives  a curve  instead  of  a 
straight  line.  Thus  the  simple  comparison  between  (I4)  and  (16) 
can  no  longer  be  made. 


The  procedure  now  outlined  was  devised  by  Lidiard 


(17) 


to 


obtain  a value  for  T without  having  to  resort  to  the  comparison 

menticned  above.  From  experimental  data  of  conductivities  measured 

as  functions  of  :i.mpu.rity  content  and  temperature  £1£2  is  plotted 

ir 


I 

—r 
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against  <T  for  various  temperatui’eso  Here  is  the  background 

impurity  content,  which  is  presei't  even  in  the  ” pure”  crystals* 

It  must  be  added  to  c to  avoid  erroneous  results  for  cases  ^ere 

c»  Cq  does  not  hold.  For  each  temperature  selected,  a smooth 

curve  is  then  dravm  through  the  plot  of  ^'*‘^0  vs*  <T  , and  from 

cr 

these  isotherms  the  values  of  cT  for  certain  selected  values  of 

c+Cq  are  taken.  Since  both  and  ^ 0“U/kT  unknown,  the  latter 

is  elimiriated  by  taking  ratios  concentratiojii  foj* 

(f  .t  jLxecTconcentration 

isotherms.  The  fixed  concentration  is  the  same  for  all  these 
ratios,  and  the  same  selected  concentrations  are  used  for  each 
isotherm.  The  only  unknown  in  these  ratios  will  be  T/T^.  We  thus 
can  obtain  sets  of  values  T/Tq  for  each  concentration,  each  value 
corresponding  to  one  temperature..  We  next  average  T/T^ofor  each 
set,  ending  up  with  a value  T/Tq  for  each  temperature  chosen.  From 
this  last  set  a mean  value  for  is  then  fourid.  Using  the  data 
of  Etzel  and  Maurer, including  a background  impurity  content 
of  Cq  = .45  X corresponding  to  10^"^  vacancies/cc,  Lidiard 

obtained  a value  of  3900*^K  for  T^.  This  corresponds  to  a binding 
energ;)r  of  .34  very  close  to  the  value  found  by  Etzel  and  Maurer. 

V/e  are  now  finally  in  a position  to  compute  the  degree  of 
association,  p,  as  a function  of  both  the  temperature  and  the 
con cnetrat ions  of  impurity  ions  used  in  the  present  experiment. 

This  is  shovm  in  Figures  6,7, Sy  and  9o  Here  p has  been  calculated 
on  the  basis  of  eq.  (7)  (dashed  lines)  and  eq.  (11)  (full  lines). 

It  is  shown  that  p falls  off  more  rapidly  with  T when  the  Coulomb 
interactions  are  taken  into  account.  This  effect  is  the  rno.^e  pro- 
noiinced  the  higher  the  impurity  content.  For  the  lowest  impurity 
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content  used  the  two  curves  of  p virtually  coincide.  It  should 
also  be  pointed  out  that  while  for  very  low  temperatures  p is 
xinity,  i,e,  all  vacancies  are  associated,  this  is  well  outside  the 
region  of  experimentation,  which  extends  from  about  250®C  to  700®C^ 
or  from  T/T^  » ,13  to  T/T^  =* 

Co  Experimental  Procedure: 

The  crystals  used  in  this  experiment  were  grown  by  means 
of  the  Kyropoulos, technique.  A platinum  crucible  of  about 
So  cm^  capacity  was  filled  with  NaCl  to  which  the  proper  amount  of 
CdCl2«2  1/2  H2O  had  been  added.  Both  reagents  were  of  the  purest 
Baker  and  Adamson  grade  available.  The  crucible  was  then  placed 
into  a one  kilowatt  open  ended  Hoskins  furnace  and  heated  slightly 
past  the  melting  point  of  NaGl  (SOO®G).  After  all  the  salt  had 
melted,  a doubla-walled  platinum  tube  with  a small  (about  3 naa  x 
3mm  X 1 cm}  NaCl  seed  crystal  attached  to  its  end  was  inserted 
into  the  melt  for  a few  seconds  and  then  vrithdravm  until  the  lower 
end  of  the  seed  crystal  was  just  beginning  to  pull  on  the  meniscus 
of  the  melt.  As  a result  of  this  dipping  procedure  a layer  of 
frozfjn  salt  was  formed  around  the  junction  of  tne  seed  crystal  and 
the  platinum  tube.  This  in  turn  provided  both  firm  mechanical 
support  and  good  thermal  contact. 

From  this  point  on  the  procedure  was  to  provide  the  proper 
thermal  gradient  to  initiate  and  continue  proper  growth  of  the 
crysta] . This  was  accomplished  by  means  of  passing  compressed  air 
through  the  platinum  tube  and  slowly  lowering  the  temperature  of 
tho  lu2-nace.  Cnee  growtdi  had  been  initiated  it  was  continued  at 
an  a\'erage  rate  of  .15  nun/min.  After  about  3 nim  were  grown  the 
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cryst-al  was  raised  until  its  3.ower  edge  was  again  pulling  on  the 
meniscus o This  procedure  was  repeated  until  a total  length  of  3 
to  4 cm  had  been  grovm.  Then  the  crystal  was  lifted  clear  out  of 
the  melt  and  allowed  to  cool  down  very  slowly  to  room  temperature. 

The  crystals  so  produced  weighed  from  about  50  to  60  grams 
and  used  up  about  half  the  melt.  It  was  also  found  empirically 
that  the  crystals  bended  to  reject  the  Cd'*"*'  ions  present  in  the 
melt.  The  ratio  of  the  Cd  added  originally  to  the  amount  found  in 
the  crystal  was  about  8 to  1.  This  ratio  was  fairly  constant 
throughout  the  entire  range  of  concentrations  used® 

The  diffusion  measurements  were  performed  in  essentially 
the  same  manner  as  described  by  >lapother,  Crooks,  and  Maurer. '"*■  * 

The  crystals  grown  as  described  above  were  annealed  and  then 
cleaved  into  specimens  about  1 cm  x 1 cm  x l/2  cm.  A piece  weighing 
about  el  gram,  cleaved  from  the  same  crystal  whose  diffusion 

coefficient  v/as  to  be  measured , was  bombarded  by  neutrons  in  the 

i 

pile  at  the  Argonne  National  Laboratory  for  about  24  hours.  During 
this  time  the  Na^^  in  the  sample  v/as  activated  to  give  (half- 

life  I4c9  hours.)  to  a sufficient  extent  to  permit  diffusion  experi- 
ments to  be  carried  on.  Upon  arrival  at  Urbana,  the  activity  of 
the  unshielded  sample  v/as  generally  found  to  be  of  the  order  of 
10  mr/hr  at  7 ft.  By  evaporation  in  high  vacuum,  a thin  layer 
(5-  6x10”^  cm)  of  radioactive  NaCl  was  made  to  form  on  one  surface 
of  each  specimen.  In  the  earlier  runs  the  specimens  were  then 

I 

' . , sealed  off,  tv/o  at  a time,  in  evacuated  pyrex  capsules.  For 

temperatures  above  550°  G the  specimens  were  first  inserted  into  a 

j 

1 stainless  steel  container,  to  pi’event  the  pyrex  from  collapsing 

t 

I 

I 

! 

! 

. ! 
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upon  thenio  Later  on  this  procedure  was  somewhat  modified,  as 
illustrated  in  Figure  2«  The  specimens  (A)  were  first  enclosed, 
again  two  at  a time,  in  platinum  capsules  (B)o  The  surfaces  bear- 
ing the  evaporated  layer  (G)  of  activated  NaCl  were  set  facing 
each  other,  kept  separate  by  a thin  platinum  spacer  (D).  The 
platinum  capsules  were  subsequently  sealea  off  in  vycor  tubes  (E) 
filled  with  helium  gas.  The  helium  for  this  purpose  was  passed 
through  a liquid  nitrogen  trap  befoi^e  use.  This  served  to  freeze 
out  most  extraneous  gases,  the  presence  of  which  might  contaminate 
the  specimens.  The  pressure  of  the  heli’^m  was  such  that  at  the 
diffusion  temperatiire  it  would  be  approximately  atmospheric.  As  a 
further  precaution  against  contamination , the  inside  walls  of  the 
vycor  tubes  had  been  coated  .vith  a thin  layer  of  NaCl  (F),  Such  a 
layer  was  prepared  by  filling  the  tubes  with  chunks  of  NaCl  and 
heating  them  to  about  600®C  under  high  vacuum,  then  cooling  to  room 
temperatureo 

From  this  point  on,  the  procedure  was  the  same  regardless  of 
whether  the  samples  were  enclosed  in  pyrex  or  vycor  capsules.  The 
capsule  was  placed  in  a stainless  steel  tube  (Figure  2,  (G)  ) 
which  could  be  closed  off  at  the  bottom.  The  temperature  could  be 
measured  by  means  of  a chromel-al.umel  thermocouple  (H).  A correction 
of  the  order  of  1®C  proved  necessary  in  the  temperatui'e  reading 
since  the  thermocouple  was  not  located  inside  the  capsule.  This 
correction  was  obtained  by  means  of  a dummy  vycor  capsule  which  was 
identical  to  the  one  shown  jji  Figure  2,  except  that  the  specimens 
were  not  radioactive  and  had  embedded,  in  them  a Pt-Pt  Rh  (10  % ) 
thermocouple o 
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For  the  actual  diffusion  runs,  the  whole  assembly  shown  in 
Figvire  2 was  placed  in  the  center  of  a fused  quartz  tube  running 
throiigh  a Hoskins  funiacea  The  temperature  of  the  furnace  was 
controlled  by  means  of  a ” Tag”  celectray  controller  to  within 
The  diffusion  times  were  adjusted  so  that  the  product  of  the  dif- 
fusion  coefficient  and  the  diffusion  time  was  approximately  lO"^  cm  .. 
However,  no  times  shorter  than  30  minutes  were  used,  because  the 
finite  rise  time  and  decay  time  of  the  temperature  would  then  be 
too  large  a fraction  of  the  total  diffusion  time.  In  order  to 
make  this  fraction  as  small  as  possible,  the  sample  was  first  put 
into  a preheating  furnace  set  at  a much  higher  temperature  for  all 
diffusion  runs  under  2 hours.  Af*.  zr  reaching  the  desired  diffusion 
temperature,  the  sample  was  then  quickly  transferred  to  the  regular 
furnace.  With  a preheating  furnace  set  at  1200°C  the  samples  could 
be  heated  from  room  temperature  to  ?00®C  within  90  seconds.  In 
addition  to  this  lower  limit  on  the  diffusion  time,  an  upper  limit 
v;as  imposed  by  the  short  half  life  of  the  radioactive  tracer  j Thus 
it  was  found  that  no  worthwhile  diffusion  measurements  could  be 
made  for  diffusion  times  exceeding  5 days,  by  which  time  the  activ- 
ity has  decayed  to  1/300  of  its  original  value. 

The  slicing  and  counting  procedure  was  identical  to  that 

(12 ) 

described  by  Mapether,  Crooks,  and  Maurer,  ' including  use  of  the 
same  microtome.  The  slices  taken  were  2 scale  divisions  (about 
4.84  X 10‘“^  cm)  thick  and  were  collected  on  scotch  tape.  Most  of 
the  slices  vrerc  counted  until  at  least  10,000  counts  had  been 
recorded,  so  that  the  error  due  to  the  distribution  of  counts  was 
held  balovi’  1%  . A number  of  corrections  haci  to  be  applied,  which 
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are  listed  as  follows: 

a)  Dead  time  of  counter  (negligible  below  4OOOO  counts/min )< 

b)  Background  (17~19)  counts/min),. 

c)  Decay  during  the  counting  of  one  slice  (negligible  for 
counts  less  than  3.5  minutes  in  duration). 

d)  Decay  during  the  counting  of  all  the  slices. 

The  corrected  number  of  counts  for  each  slice  in  a specimen 
was  then  plotted  semilogarithmicaliy  against  the  square  of  the 
depth  X.  Unfortunately,  the  advance  of  the  knife  assembly  down- 
ward was  not  urJjform,  so  that  x could  not  be  taken  as  n times  the 
thickness  of  one  slice,  where  n is  an  integer.  Through  calibration 
with  the  interferometer  mounted  in  the  knife  assembly  of  the  micro- 
tome the  exact  depth  at  each  slice  was  determined.  Thus  x was 
still  taken  as  n times  the  tliickness  of  one  slice  (set  arbitrarily 

I 

at  4.84  X 10"‘*  cm  for  all  slices)  but  n now  was  no  longer  a set  of 
integers,  but  depended  on  the  vertical  position  of  the  knife.  A 
typical  ■''emilogarithmic  plot  of  counts  I vs . is  shown  in  Figure  .5. 

From  the  slope  of  such  a plot  as  well  as  a knowledge  of  the 
diffusion  time  the  diffusion  coefficient  can  be  determined.  The 
equation  (derived  in  appendi.c  II)  is  as  follows 

0 ^ .X.,  (19) 

4t  Inl^-lnla 

vrfiere 

t = time  of  diffusion  in  seconds 

= activity  in  counts/min  at  point  a 
' activity  in  counts/min  at  point  b 
n^  = depth  in  crystal  at  point  a in  units  of  one  slice 
thlclmess 
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* depth  in  crystal  at  point  b in  units  of  one  slice 
thickness 

\ = thickness  of  one  slice  in  cm 
D = diffusion  coefficient  in  cm^/ 

The  errors  involved  in  the  diffusion  coefficient  arise 
from  the  following  sources; 

a)  Errors  in  diffusion  time.  These  errors  are  appreciable 
only  for  very  short  times  and  negligible  for  times  longer  than  a 
few  hours. 


b)  Errors  in  the  slope,  i,e,  , In  most  cases. 

Inlb-lnl^ 

these  errors  are  only  a feiv  percent.  However,  for  diffusion 
coefficients  in  one  range  of  10“^^  to  lO"^^  cmVsec,  where  only  a 
few  countable  slices  can  be  taken,  these  errors  become  very  large. 

c)  Errors  in  depth  measurement.  Such  errors  are  due  to 
the  fact  that  nX  does  not  represent  exactly  the  depth  of  a slice. 
This  is  because  of  errors  in  the  calibration  and  a slight  change 


in  calibration  of  the  downward  movement  of  the  knife  assembly  when 
it  moves  through  the  same  vertical  range  again.  These  errors  are 
estimated  to  be  about  2%  , but  will  contribute  a 4 9©  error  to  D, 
because  the  depth  appears  as  the  square  in  eq.(19). 

d)  Errors  arising  from  approximations  made  involving 
solution  of  tne  diffusion  equation.  This  is  discussed  more  fully 
in  appendix  II.  At  this  point  it  may  be  stated,  hov.’svsr.  that 
these  errors  are  mainly  due  to  not  using  the  exact  solution  of  the 
diffusion  equation  and  not  taking  into  account  the  finite  thickness 

sjx  ^ fciXwTo  dx  c:  uulUtJ  »Ulci-LX  X Or  High  QllfuSiOn 

coefficients  but  increase  very  rap>idly  for  low  oneSu 
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e)  Errors  in  temperature.  It  is  believed  that  the  tempera- 
txire  is  known  to  within  l®Co  This  causes  an  error  in  D ranging 
from  about  1 in  the  high  temperature  region  to  4%  at  the 
lowest  temperatures  used  in  this  experiment. 

f)  All  other  errors.  They  are  believed  not  to  exceed  1 /q» 
Since  some  of  the  errors  depend  directly  or  indirectly  on 

the  magnitude  of  the  diffusion  coefficient,  we  will  list  these 
errors  in  table  I for  various  values  of  D. 


Table  I 


D(cm^/sec)  Source  of  error  in  D(  ^ ) Total  error 

ln  b(  Ijr 

Time  Slope  Depth  Soln.of  Temp.  All  other 


equat ion 

10"^  5 2 4 3 1 1 16 


10-9 

2 

2 

4 

10-10 

3. 

2 

4 

10-11 

- 

3 

4 

10"!^ 

- 

15 

4 

10-13 

30 

4 

3 

1 

1 

13 

3 

1 

1 

12 

7 

2 

1 

17 

12 

3 

1 

35 

25 

4 

1 

64 

In  general,  then,  we  can  say  that  the  errors  in  the  measured 
diffusion  coefficient  are  about  15  % for  all  but  very  low  values  of 
D,  and  that  for  these  values  the  accuracy  of  D diminishes  very  fast. 
The  conductivity  was  measured  by  a method  similar  to  that 

/ T T \ 

used  by  Bean.^  ' A sample  was  prepared  by  first  cleaving  a plate 
about  1 mm  phick  from  a diffusion  specimen  prior  to  diffusion. 

Then  both  sides  were  painted  with  a suspension  of  graphite  in 
alcohol,  after  v/hich  the  actual  sample  (about  2 mm  x 2mm)  was 


f*' 

I 


i 
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cleaved  outo  The  sample  was  t-hen  placed  between  the  two  electrodes 
or  a crystal  holder o The  original  holder  used  in  this  experiment 
was  identical  to  the  one  described  by  Bean,  its  main  feature  con- 
sisting of  two  1 mil  platinum  electrodes  spot-v/elded  onto  n5.ckel 
buttons  at  the  end  of  a thin-walled  stainless  steel  tubeo  During 
the  course  of  the  experiment  it  was  decided  to  modify  this  slightly 
to  avoid  possible  contamination  of  the  crystal  due  to  the  presence 
of  nickel  and  stainless  steel  at  high  temperatures.  The  improved 
holder  is  shown  in  Figure  3.  It  consists  of  two  solid  platinum 
electrodes  (A),  one  attached  to  the  end  of  vycor  tube  (B),  the  other 
one  moving  freely  inside  another  vycor  tube  (Cj,  but  held  in  place 
by  3 platinum  wires  (only  2 of  which  are  shown  in  the  Figure)  and 
springs  (D)  in  the  cold  temperature  region.  A Ft-Pt  Rh(10%  ) 
thermocouple  (E)  is  led  in  throi2gh  two  small  Stupakoff  seals.  The 
platinum  wires  holding  the  lower  electrode  constitute  the  ground 


lead.  The  high  lead  consists  of  a platinum  wire  (F)  running  inside 
the  inner  vycor  tube  and  being  led  out  of  the  assembly  by  means  of 
another  Stupakoff  seal  (C)  and  an  amphenol  socket  (H). 

After  a sample  (I)  had  been  put  into  place,  the  whole 
assembly  was  inserted  into  a vacuum  tight  quartz  tube  running 
through  a Hoskins  furnace  and  helium  continuously  circulated  through 
the  tube.  The  furnace  temperature  was  regulated  by  a " Tag*^ 
celectray  controller  to  within  t Care  was  taken  during  each 

conductivity  rim  to  adhere  as  closely  as  possible  to  the  same 
thermal  cycle.  This  consisted  cf  taking  readings  about  25  minutes 


noart  at  SOPCifipH  t.  Pmnor’.'at-.IIV'O  -iri-t- ® 


or  AOr>  ^ j irtr\r\Or% 
\J  duu  /UU  O. 


the  keeping  the  sample  at  700*^’c  for  about  90  minutes,  and  then 
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taking  readings  at  the  same  time  and  teraperatinre  intervals  while 
the  temperature  was  lowered  again® 

To  measure  the  conductivity,  an  AoC,,  bridge  technioue  was 
used.  The  circuit  is  illusti*ated  schtinat ically  in  Figure  4«  The 
sample  holder  (C)  is  connected  i.ito  one  arm  of  a Gexiaral  Radio 
716-C  bridge.  The  generator  is  a Hewlett-Packard  audio-oscillator 
(A),  the  null-meter  a General  Radio  736-A  v;ave  analyser  (D).  In 
addition,  a General  Radio  721-D  precision  condenser  (E)  has  been 
put  into  the  circuit  for  balancing  purposes.  Resistances  from 
10-“^  ohms  doivn  to  lO'^  ohms  were  measured  by  a substitution  method, 
balancing  Gq  and  Cj^  with  and  v/ithout  the  sample  in  the  circuit. 

The  frequencies  used  ranged  from  100  cycles/sec  for  the  highest 
resistances  to  10000  cyclea/^ec  for  the  lowest  ones,  most  of  the 
readings  being  taken  either  £.t  1 kc  or  3.0  kc.  The  accuracy  of  the 
resistance  measurements  is  about  4%,  including  a 1^  error  in 
the  frequency,  for  resistances  iov/er  than  10^  ohms,  an  eoual 
bridge  arm  technique  v/as  usee.  In  this  arrangement  the  upper  arms 
of  the  bridge  were  adjusted  to  be  equal  and  a Shallcross  px'*ecision 
resistance  box  placed  in  parcllei  with  the  external  balancing 
condenser.  The  value  of  the  resistance  needed  for  balance  gives 
directly  the  parallel  resistance  of  the  sample  and  its  holder. 

The  resistance  box  was  calibrated  and  found  to  be  accurate  to  better 
than  1 

In  computing  the  total  errors  in  the  conductivity,  we  must 
take  Into  account  two  other  sources  of  error  in  addition  to  the 
ones  nentioried  above.  One  of  thF.se  is  the  uncertainty  in  the 
tenipei'atur e , which  causes  exTors  in  the  conductivity  ranging  from 
4%  at  250^3  to  1 at  700*^0,  ‘‘he  other  source  is  a slight  change 


in  the  dine  nions  of  the  ;an.ple  at  high  temperatures  due  to 
evaporation.  This  contrit.ites  at  most  a 2 % error,  ana  can  be 
neglected  at  low  temperatuves.  Vtfe  thus  can  state  that  the  total 
error  in  the  conductivity  r nges  from  about  6 ^ at  low  temperatures 
(high  resistances)  down  to  / % at  high  temperatures.  The  con- 
ductivities are  therefore  kr  oim  wLth  much  greater  precision  than 
the  diffusion  coeff icientf}. 
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Do  Results: 

Foia*  cryjjials  of  NaC.  with  varying  concentrations  of  CdCl2 

were  used  in  the  conductivity  and  diffusion  measurements.  Various 

sections  of  these  crystals  vere  analyzed  polarographically  for 

Cd  content.  In  table  II,  ths  minimum,  maximum  and  average  molar 
Ccl 

ratios  are  listed  for  theie  crystals.  It  should  be  po:jited  out 
that  since  some  sections  of  tie  crystals  were  used  more  tl.an  once, 
both  for  conductivity  and  difJusion  measurements,  the  avei’age  is 
weighted  to  reflect  not  only  tie  concentration  of  each  section  but 
also  the  number  of  times  it  wa.  used. 


Tab.e  II 


crystal 

minimum  ^ x 10^ 

riixinmm  ^ x 10^ 
TTa 

average 

A 

26 

33 

3<» 

B 

11 

1^ 

It 

C 

3 

e 

✓ 

1. 

D 

t-4 

o5 

The  accuracy  of  the  polaro.vraphic  analysis  ranges  from  about  25®/^ 
for  the  lov/est  concentration  to  5 °o  highest. 

Figure  1C  shows  the  result  of  the  conductivity  mec.surements 
on  the  four  crystals  contaojiing  Cd  I curves  A,  B,  C,  and  D , and 
also,  for  comparison,  the  <;onductivity  of  a typical  ” purvi**  NaCl 
crystal  (curve  S),  These  data  in  eac)i  case  represent  an  average 
taken  from  4 to  6 different  runs,  some  of  which  were  taker  on 
different  sections  of  the  crystal  and  some  repeated  on  the  same 
section.  Each  run  was  taken  through  the  same  thermal  cyc]e,  as 
outlined  iii  the  preceding  sectiono  A significant  decreast  has  been 
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observed  in  the  conductivity  of  crystals  A,  B,  C,  and  D during  the 
second  half  of  the  cycle.  This  decrease  is  the  more  pronounced  the 
higher  the  concentration,  ranging  from  about  for  the  lowest  to 

40  for  the  highest.  No  such  decrease  has  been  observed  for 
**pure”  NaCl  crystals  or  crystals  containing  instead  of  Cd. 

Similar  observations  have  also  been  made  by  Btzel  and  Insurer. 

The  evidence  based  on  their  as  well  as  the  recent  work  Is  far  from 
conclusive.  However,  it  appears  that  both  a volume  and  ai  elec- 
trode-crystal interface  effect  are  present.  In  the  forme:.*,  Cd'*"*’ 
ions  are  somehow  removed  from  the  volume  of  the  cx*ystal,  >erhe.ps 
to  form  larger  ag;gregates.  This  would  effectively  cut  dom  the 
number  of  vacancies  and  hei^ice  the  conductivity.  As  far  as  the 
second  effect  is  concerned,  it  is  believed  that  some  form  of  a high** 
resistance  layer  is  formed  at  the  interface  between  the  crystal  and 
the  graphite  electrode.  It-  is  not  clear  why  this  should  occur  when 
Cd'*"*'  ions  are  present  and  not  when  we  are  dealing  with  tho  same 
experimental  arrangement  but  with  Ca'*"*'  Ions  present  or  eliie  no 
deliberately  added  impurity.  The  whole  problem  needs  to  be  investi- 
gated in  greater  detail.  Xt  can  be  stated,  however,  that  on  the 
basis  of  the  present  evidence  it  appears  that  the  volume  effect 
outweifiis  any  other  effect  c 

Figures  11,  12,  13,  and  I4  show  the  diffusion  coefficients 
measured  by  means  of  radioactive  tracer  methods.  On  the  same  plot 
the  diffusion  coefficients  calculated  from  the  measured  conduc- 
tivities by  means  of  the  Einstein  relation  are  shown  (full  lines). 
The  conductivities  used  for  this  purpose  are  the  ones  shovn  in  the 
lower  half  of  each  cycle  in  Figure  10,  \Ie  have  also  calculated 
the  diffusion  coe.fficients  based  on  what  the  ccwiductivitics  would 
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be  were  there  no  association.  The  relation  used  was 

Measured  * ^no  association  (20) 

This  is  also  shown  In  Figures  11,  12,  1?.  and  14»  The  dotted  lines 
correspond  to  p calculated  from  (7),  and  the  dashed  lines  to  p 
calculated  from  (11).  In  Figure  14  these  two  lines  are  so  close 
together  as  to  virtually  coincide. 

As  expected,  the  curves  based  on  lie 

no  association 

higher  than  the  ones  based  on  Further,  since  p 

measured  * 

determined  from  (7)  is  greater  than  p determined  from  (11),  the 
values  of  0 calculated  based  on  the  simple  theory  of  association  are 
larger  tha;i  those  based  on  the  more  detailed  theory.  It  can  also 
be  noted  from  these  Figures  that  the  effect  of  taking  p into  account 
increases  with  concentration  and  decreases  with  temperature. 

The  failure  or  satisfaction  of  the  Einstein  relation  and  the 
conclusions  that  can  be  drawn  from  it  will  be  discussed  in  the  next 
section.  V/e  wdll  limit  ourselves  at  this  point  to  some  remarks  an 
the  scattering  of  the  measured  diffusion  coefficients.  As  already 
pointed  out  in  the  preceding  section,  the  diffusion  coefficients 
obtained  by  radioactive  tracer  techniques  are  subject  to  much 
greater  errors  than  the  conductivity  measurements.  In  addition  to 
these  errors,  the  spread  in  concentration  of  impurity  over  a crystal, 
in  view  of  the  relatively  large  size  of  the  diffusion  specimens, 
causes  errors  of  up  to  20  % , which  is  of  about  the  same  oirder  as  all 
the  other  errors  combined.  An  attempt  has  been  made  to  minimize 
these  additional  errors  by  cleaving  the  specimens  from  as  small  a 
section  of  the  cx'ystai  as  possible,  and  by  using  them  over  again 
two  or  three  times  at  different  temperatures.  Nevertheless,  in 
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view  of  the  large  errors  and  the  small  number  of  points  involved, 
as  compared  to  conductivity  measurements,  an  averaging  process  is 
probably  meaningless  and  has  been  omitted  here*  The  scattering  of 
the  measured  diffusion  coefficients  shown  in  Figures  11,  12,  13, 
and  14  is,  \irith  one  or  two  exceptions,  well  within  the  range  of 
experimental  error*  IVhere  measurements  on  two  specimens  at  the  same 
temperature  were  made,  the  larger  diffusion  coefficients  usually 
corresponds  to  a specimen  with  a slightly  higher  Cd  content*  But 
even  here,  in  a few  instances,  the  order  is  reversed,  indicating 
that  the  other  errors,  discussed  in  the  previous  section,  are  also 
appreciable.  It  should  be  pointed  out,  however,  that  there  is  less 
scattering  of  diffusion  data  in  the  present  work  than  has  been  found 
in  the  intrinsic  region  of  ’’pure**  salts,  such  as  NaCl,^^^*^ 

and  even  though  presumably  there  were  no  concentration 

gradients  in  these  "pure**  ciystals* 

T?  n O G O 'I  r\T\ 

\ie  first  consider  the  conductivities  shown  in  Figure  10* 

In  the  low  temperature  region  (below  about  400®C)  they  fellow  more 

or  less  strai^t  and  parallel  lines.  The  slope  is  given  by  -U/kT^, 

where  U is  the  activation  energy  for  the  jump  of  a vacancy.  The 

value  of  U obtaiiied  from  the  present  data  is  ,&3  eV,  the  same  /.  s 

(10) 

was  found  by  Etzel  and  Maurer* ' 

This  very  uniform  behavior  of  the  conductivities  in  the  low 
temperature  region  caai  easily  be  explained  by  the  fact  that  we  are 
entirely  in  the  stmicture  sensitive  region.  This  means  that  almost 
all  of  the  vacancies  contr’ibuting  to  the  conductivitj^  are  those 

introducce:  by  jjmpurities.  Hence  the  energy  of  creation  of  a vacancy 
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does  not  enter  into  the  expression  for  the  conductivity.  Strictly 
speaking,  there  is  a temperature  dependence  in  the  pre-exponential 
factor  of  the  conductivity,  i.e, 

<r-  * e -“/MJ  (21) 

However,  the  change  in  \-/ith  T in  the  low  tempei-atui*e  region  is 
aa.most  entirely  due  to  the  exponent. 

In  the  high  temperatui'e  region  the  conductivities  curve 
upward,  as  expected,  since  now  the  thermally  generated  vacancies 
begin  to  equal  and  eventually  outnumber  the  impurity  vacancies^ 

Thus  we  can  write 

(22) 

where  V/  is  the  energy  necessary  to  create  a vacancy  pair.  It  is 
further  apparent  from  FigtjLre  10  that  the  lower  the  icipurity  con- 
centration, the  lower  the  temperature  at  which  the  upward  curving 
of  the  conductivities  begins. 

This  simple  treatmait,  however,  fails  to  explain  why  the 
slopes  are  not  constant  until  the  thermal  vacancies  begin  to  con- 
tribute. At  temperatures  well  below  the  ones  where  thi»5  contri- 
bution appears  tliere  is  a distinct  downward  curving  of  the  conduc- 
tivities, as  can  be  seen  from  Figure  10,  A possible  explanation 
of  this  phenomenon  has  been  suggested  by  .According  to 

his  theory,  the  jump  frequency  of  a vacancy  will  be  affected  by  the 
vibrations  of  the  two  negative  ions  that  stand  in  the  way,  since  it 
is  easier  for  the  vacancy  to  jump  when  the  two  negative  ions  are 
moving  away  from  each  other.  Thus  the  pre-exponential  factor  in 
(21)  is  changed.  According  to  Bean,  the  result  is 


1 

I 

I 

I 

I 

1 
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<r.  K (23) 

T-’ 

The  present  data  are  indeed  fitted  better  by  (23)  than  by 
(21),  but  the  fit  is  by  no  means  perfect,  A good  fit  would  result 
from  in  the  denominator  of  (23)  instead  of  T^,  Thus  other 

factors  must  be  operating  in  the  system  NaCl  4-  CdCl«  for  which  at 

. ! 

present  there  is  no  simple  explanation. 

Before  passing  to  a comparison  of  the  measured  and  calcu- 
lated diffusion  coefficients  shown  in  Figures  11,  12,  13,  and  14j 
let  us  first  consider  some  of  the  assumptions  involved.  The  first 
assumption  made  is  that  the  same  volume  effect  which  causes  a 
decrease  in  the  conductivity  also  causes  a decrease  in  the  diffusion. 

Since  both  are  based  on  a vacancy  mechanism  and  since  the  net  result 

of  the  voliime  effect  is  to  cut  down  the  number  of  vacancies,  such 

an  assumption  se  ciuw  x^ccscrx  able.  It  is  for  this  reason  that  in  com-  | 

paring  the  measui'sd  and  calculated  diffusion  coefficients  the  lov:sr 

half  of  the  conductivity  in  each  cycle  was  used.  The  probable 

correctness  of  this  assumption  is  further  emphasized  by  the  fact 

that  using  it  ths  Einstein  relation  is  satisfied  In  the  hi^ 

temperature  region,  Vvhereas  using  the  higher  values  of  d’  (upper 

half  of  cycle would  lie  well  above 

there  have  been  a number  of  mechanisms  proposed  vrtiich  would  cause 

^calculated  ^measured » permit  the 

opposite.  Indeed,  previous  work  on  self-diffusion  in  " pure”  alkali 
halides  (12,13)  j^^s  always  either  resulted  in  the  Einstejji  relation 
being  satisfied  or  else  some  mechanism  has  been  proposed  which  would 
permit  some  vacaiicies  to  contribute  to  diffusion  but  not  to  conduc“ 
tiv3  ty. 

■ i 

I 
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The  second  assumption  concerns  the  effect  of  the  atmosphere 
of  impurity  ions  in  hiiidering  the  motion  of  vacancies  (see  eq.  (l8)  )• 
If  one  fixes  a point  of  reference  inside  a crystal,  most  of  the 
motion  of  the  vacancies  past  this  point  under  the  influence  of  a 
field  v/ould  appear  to  be  a random  jumping  rather  than  drift.  In 

the  steady  state,  however,  a definite  displacement  of  the  vacancies 
with  respect  to  the  impurity  ions  will  have  occurred  under  the 
influence  of  a field,  v/hile  it  would  not  be  reasonable  to  assume 
that  such  a displacement  occurs  when  no  field  is  present.  Thus  the 
retardation  effec:  is  present  in  the  conductivity  and  not  in  the 
diffusion.  Therefor"  in  applying  the  Einstein  relation  ^f^^cu^^ated 
should  be  corrected  by  the  factor  (1-f ) while  ^measured  sl-^uld  not. 
This  has  been  neglected  here,  because  the  correction  even  for  the 
hipest  impurity  concentration  would  only  be  a few  percent,  and 
quite  negligible  for  lower  concentrations.  For  a much  hig;her 
concentration  of  impurity,  on  the  other  hand,  such  as  can  be  achieved 
with  Ca  and  perhaps  other  metals,  though  not  Cd,  this  effect  could 
not  be  neglected. 

The  third  assumption  concerns  the  binding  energy  of  a com- 
plex, v.'hich  v.»e  have  taken  to  be  ,34  eV,  This  value  was  aiTived  at 
by  fixing  both  the  binding  energy  of  a complex  and  the  ma^^nitude 
of  the  interactions  at  large  distances.  This  is  obviously  a much 
better  value  than  the  one  obtained  from  eq,  (12),  i,e,  ,60  eV, 
completely  neglecting  these  long-range  interactions.  Furthermore, 
using  Tp  " 7000®K  in  eq,  (7)  or  eq,  (11)  would  give  degretjs  of 
association  at  or  near  j.OO  for  most  of  the  crystals  used  in  the 
present  work,  an  obviously  absurd  conclusion „ However,  by  taking 
o34  sV  as  the  reel  bindjng  energy  of  a complex,  we  underestimate 


somewhat  the  effect  of  the  long-range  interactions.  This  in  turn 
results  in  overestimating  the  degree  of  association,  calculated  by 
eq.  (11),  In  order  to  apply  the  theory  correctly  it  is  necessary 
to  take  into  account  the  non-Coulombic  Interactions  between  the 
vacancies  and  the  impurity  ions,  both  at  nearest  neighbor  sepa- 
rations and  at  3.arger  distances.  These  interactions  are  due  to  the 
repulsive  energj-,  distortion  of  the  lattice  around  a defect,  and 
polarization  due  to  dipoles  created  by  displacement  of  ions  and 
electronic  polai-izability . They  have  been  considered  previously 
by  Bom  and  Meyer,  ^21)  j^ott  and  Littleton,  (7)  and  Reitz  and 
Gammel.^22)  letter  obtained  a binding  energy  for  the  Cd'^^  ion- 

vacancy  complex  in  NaCl  of  .44  eV,  vdiich  is  almost  the  same  as  the 
Coulomb  interaction  energy  of  an  impurity  ion  and  a vacancy  in  next 
nearest  neighbo.r  positions  (.43  eV),  A very  recent  computation  by 
Bassani  and  contains  some  improvemoDts  over  the  work  of 

Reitz  and  Gammel.  Their  result  for  the  binding  energ>-  of  th«  complex 

in  ques1;ion  is  Vno  Tr  V\  A £%  ^ ^ l_  _ A ^ ^ 

uu  uiie  uae  present 

theoretical  value,  Bassani  and  Fumi  did  consider  non-Coulombic 
intei'actions  as  far  rway  as  next  nearest  neighbor  positions.  Now 
Lidiard^^7)  points  out  that  it  does  not  make  much  difference  in 
connection  with  our  problem  whether  we  take  the  limiting  distance 
differentiating  between  complexes  and  unassociated  charges  to  be 
twice  the  lattice  constant  or  ten  times  the  lattice  constant.  If 
we  can  therefore  assume  that  non-Goulombic  terms  are  of  little 
importance  past  the  next-nearest  neighbor  separation,  to  apply  the 
theory  correctly  we  need  only  use  ,36  eV  instead  of  ,34  oV  for  the 
binding  energy  and  assume  the  limiting  distance  of  association  to 
be  the  next-nearest  neighbor  distance.  This  would  result.  In  p 


( 
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Increasing  slightly  at  any  given  temperature,  regardless  of  whether 
eq*  (7)  or  eq«  (11)  is  used.  The  only  correctioi  still  necessary 
to  the  theory  of  association  would  now  be  the  non-Coulomb ic  inter- 
actions at  larger  than  next-nearest  neighbor  positions.  These  are 
presxuned  to  be  quite  small. 

The  main  result  of  our  assumption  regarding  the  binding 
energy  is  then  that  a small  correction  may  be  necessary  vrhich  would 
shift  the  curves  of  on  either  eq.  (7)  or  eq.  (11) 

sli^tly  upward.  It  would  not,  however,  materially  affect  any  of 
the  conclusions  reached  as  a result  of  this  work. 

We  are  now  in  a position  to  discuss  the  main  part  of  the 
present  work,  that  is,  the  failure  or  satisfaction  of  the  Einstein 
relation  as  a function  of  concentration  and  temperature  and  the 
con clxis ions  that  can  be  drawn  there  from  regarding  the  degree  of 
association.  From  the  results  shown  in  Figures  11,  12,  13,  and  14 
it  is  apparent  that  the  Einstein  relation  is  satisfied  for  all  con- 
centrations in  the  high  temperature  region.  However,  since  the 
degree  of  association  even  for  the  highest  concentration  is  <xiiy 
about  10%  at  these  temperatures,  it  is  impossible  to  say,  in  view 
of  the  large  eiror  In  , ^diether  the  Einstein  relation  holds 

on  the  basis  of  Pleasured  ^no  association*  ^ 
temperature  region,  however,  both  crystals  A and  B show  a signifi- 
cant deviation  from  the  Einstein  relation  if  we  use  ^measured 
a basis.  This  deviation  is  well  outside  the  limit  of  experimental 
error.  This  is  believed  to  be  definite  evidence  that  there  is  a 
considerable  degree  of  association  in  these  crystals  in  the  low 
temperature  region.  iVhen  we  compare  “calculated 

on  the  basis  of  eq.  (7)  and  eq.  (11),  we  find  that  for  crystal  a 
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•qr  (7)  represents  the  degree  of  association  some^rtiat  better  than 
eqo  {11).  In  crystal  B the  evidence  is  anich  less  conclusive, 
especially  since  the  curves  representing  basis 

of  the  two  equations  lie  closer  together.  There  the  errors  in 
^measured  same  order  as  the  difference  between  these  two 

curves.  However  there  is  no  question  that  a sigjiificant  degree  of 
association  also  is  present  in  crystal  B. 

In  crystals  C and  D no  deviation  from  the  EJnstein  relation 


can  be  found.  But  here  on  the  basis  of  o^measured 

either  form  of  association  within  the  range  of  probable 

error  of  so  this  behavior  is  understandable.  There  is, 

however,  a significant  point  that  must  be  brought  up  here  in  this 

(12 ) 

connection.  The  work  of  Mapother,  Crooks,  and  Maurer'  has  shown 
that  a significant  deviation  from  the  Einstein  relation  occurs  in 
♦*pure”  NaCl  and  NaBr  in  the  low  temperature  reglcn*  This  deviation 
correpsonds  to  a degree  of  association  of  about  50^/p*  It  is  believed 
that  this  very  large  degree  of  association  is  due  to  one  or  more 
types  of  impurity  ions  in  the  background  impurity  content,  vdiieh, 
in  contrast  to  Cd,  are  very  highly  associated,  Spectrographic 
analysis  of  »*pure"  NaCl  shoves  at  least  2j  impurity  elements  present. 
Of  these  Fe,  Mn,  Mg,  As,  Sn,  Sb,  Hg,  Pb,  Cu,  Zn,  Mi,  Co,  Ba,  Ca, 

Cr,  Al,  Zr,  and  Ti  are  possible  candidates  for  such  an  effect.  Let 
us  now  on  the  other  hand  consider  crystal  0.  The  measured  conduc* 
tivity  is  about  six  times  and  the  measured  diffusion  coefficient 
three  times  that  of  liapother’s  ” pure”  NaCl  crystals.  In  a crystal 
containing  2N^  Cd  ions  wa  would  thus  have  Nj^  of  these  highly 
associated  background  impurity  ions  to  obtain  the  same  ratios  for 


the  diffusion  coefficients.  But  while  the  2N^^  Cd  atoms  would 
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contribute  equally  to  the  conductivity  (the  degree  of  association 

being  negligible)  only  half  the  background  impurity  ions  would 

contribute,  assurr.ing  a 50  % degree  of  association.  Thus  for  crystal  D 

we  would  expect  j^measured  ■ 3 ■ 1 » Hence  p *1/6  or 

® calculated  2 1/2  1-P 

16.6  yo  • This  would  certainly  not  cause  a deviation  beyond  the 
limit  of  error  in  ^measured  very  low  diffusion  coefficients 

In  crystal  D,  and  therefore  it  is  not  surprising  that  the  effect  of 
the  background  iiapurities  cannot  be  detected.  This  illr’stration, 
hovever,  serves  to  show  that  for  higher  concentrations  of  Cd  the 
eff  ect  of  the  background  impurities  will  be  completely  masked. 

In  summary,  then,  we  may  say  that  a deviation  from  the 
Einstein  relation  can  be  detected  In  those  crystals  where  it  is 
large  ei  ough  to  be  beyond  the  range  of  the  experimental  error  in 
^measured*  This  deviation  is  approxiuately  that  expected  on  the 
basis  of  the  theory  of  association.  However,  it  is  not  possible 
to  state  coacliasively  whether  the  present  data  are  fitted  better  by 
the  simple  theory  of  association  (eq.  (7)  ) or  by  the  more  detailed 
theorj--  resulting  in  eq.  (11). 

In  cona'  usion,  it  might  be  of  value  to  suggest  a few  related 
experiments  \-rtiich  might  shed  further  light  on  the  subject  of 

association.  This  would  also  include  improving  the  present  experi- 
ment by  producing  a more  uniform  Cd^*  ion  concentration  (higher 

ratio)  and  if  possible  also  crystals  of  hi^er  Cd  content, 

crystal 

The  former  would  cut  down  the  error  in  Dmeasvir ed * the  latter 

would  produce  a greatei-  degree  of  association,  which  mi^t  be 
checked  more  quantitatively.  Next,  a different  divalent  ion,  like 
Ca’’"*'  or  St*'*' g and  a different  solvent  salt,  like  KCl,  might  be  used. 
Crystals  of  a molar  ratio  Ca/Na  up  to  25  x 10“^  have  been  grown  by 
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(11)  (23) 

Bean*  Since  accordLng  to  Bassani  and  Funsi  the  Ga,  Gd, 

and  Sr  complexes  in  both  NaCl  and  KCl  have  approximately  the  same 

binding  energy  (all  lie  between  ,32  eV  and  ,45  eV),  we  would  expect 

very  large  degrees  of  association  for  high  impurity  content.  Thus, 

it  may  be  possible  to  decide  quantitatively  whether  eq,  (?)  or  eq, 

(11)  better  describes  the  degree  of  association.  In  addition,  the 

system  NaCl  ♦ CaCl2  has  the  advantage  of  not  having  the  annoying 

decrease  of  conductivity  with  time.  Going  now  to  pure  crystals, 

the  theory  that  highly  associated  background  impurity  ions  are 

(12 ) 

responsible  for  the  results  of  liapother.  Crooks,  and  I-laurer, 
could  be  checked  by  producing  *»pure**  crystals  (possibly  by  zone 
melting)  that  have  one  or  two  orders  of  magiitude  less  Impurity 
content  than  the  usual  **pure*’  crystals  (10^*^  vacancies/co).  There 
should  be  much  less  association  in  such  very  pure  crystals.  The 
disadvantage  of  having  very  low  diffusion  coefficients  with  a con- 
sequently large  probable  error  may  make  this  experiment  difficult, 
however- 

Finally,  the  diffusion  coefficient  of  a divalent  impurity 
ion  may  be  measured  directly  by  using  an  isotope  of  this  ion  as  a 
tracer.  If  a large  degree  of  association  exists,  then  such  an  ion, 
having  so  to  speak  a built-in  vacancy  into  which  it  can  jiimp,  should 
have  an  exceedingly  high  diffusion  coefficient.  The  difficulty 
here  lies  in  activating  the  relatively  few  impurity  ions  In  the 
lattice  to  a sufficient  extent  to  carry  out  the  tracer  technique. 
Such  an  experiment  has  failed  for  this  reason  with  the  system 
NaCl  -V  GuGl2o  It  may  succeed  with  the  system  NaCl  — CaGl.^,  since 
10  times  more  Ca  than  Gd  can  be  bviilt  substitutionally  into  the 
la'^tice.  Also,  t.lii  slow-nautron  cross-section  of  Ca  is  considerably 
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hi^er  than  that  of  Cd  for  the  isotope  suitable  for  a diffusion 
experiment.  Success  of  such  an  experiment  mi^t  be  insured  by 
making  use  of  a neutron  flux  higher  than  is  now  available  either 
at  Oak  Ridge  or  at  the  Argonne  National  Laboratory. 
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CHAPTER  II 

SELF  DIFFUSION  IN  PURE  POTASSIUM  CHLORIDE 


A.  Introduction; 

While  the  work  described  in  Cht.pt er  I was  in  progress, 

H,  Witt  of  the  University  of  Goettingen  published  data^^^^^  on  the 

temperature  dependence  of  the  self  diffusion  coefficient  of  the 

potassium  ion  in  potassium  chlorideo  In  the  neighborhood  of  600®C, 

^ere  the  vacancy  content  of  the  crystals  should  be  uninfluenced 

by  impurities,  Witt  found  a deviation  from  the  Einstein  relation 

between  diffusion  coefficient  and  electrical  conductivity.  This 

behavior  of  KCl  is  unlike  that  of  NaCl  and  NaBr  as  determined  by 

(12  ' 

the  work  of  i”apother,  Crooks,  and  i-iaurer,'  * As  already  indicated 
in  the  previous  chapter,  at  temperatures  sufficiently  high  that  the 
accidental  impurities  had  no  effect  upon  the  vacancy  concentration, 
these  investigators  found  the  Einstein  relation  to  be  satisfied  in 
NaCl  and  NaBr, 


The  present  work  ’.-.'■as  undertaken  to  check  the  data  of  Witt 
since  his  results  had  rather  interesting  implications.  In  par- 
ticular, Witt’s  data  furnish  evidence  for  the  existence  of  pairs 
of  positive  and  negative  ion  vacancies  in  KOI.  Such  pairs  are 
illustrated  in  Figure  IF,  The  mechanism  by  which  pairs  may  account 
for  failure  of  the  Einstein  relation  has  already  been  discussed 
(see  p,3)o 


•There  is  some  additional  evidence  suggesting  although  not 
conclusively  establishing  the  existence  of  pairs,  Seitz has 
proposed  that  pairs  may  be  responsible  for  the  relatively  rapid 
migr'ation  of  color  centers  in  KGl,  as  observed  by  Estermann.^  Leivo., 
and  Sbern,  and  for  the  enhanced  darkening  of  X-ray  colored  KCl 
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/ 27  28 ) 

and  NaCl  with  additive  CaGl„,  ’ Also,  the  formation  of  pairs 

(25 ) 

is  favored  by  the  fact,  pointed  out  by  Seitss,'  that  there  is  no 
energy  miniiaum  for  the  configuration  of  vacancies  of  opposite  Jign 
in  next-nearest  neighbor  positions.  This  means  that  such  vacancies 
■will  automatically  be  drawn  together  to  form  a pair  when  they 
approach  as  closely  as  the  next-nearest  neighbor  distance. 

In  viev/  of  all  these  facts  it  was  therefore  deemed  advis- 
able to  determine  as  accurately  as  possible  v/hether  pairs  contribute 
appreciably  to  the  self  diffusion  coefficient  of  the  cation  in  KCl. 
VvTiile  a negative  answer  would  not  nile  out  the  existence  of  pairs, 
an  affirmative  result  would  sake  more  'plausible  the  explanation  of 
other  phenomena  in  nerms  of  the  pair  hypothesis. 


Be  Experimental  Pi-ocedure: 

With  a few  modifications,  the  experimental  procedure  for 
KCl  vras  the  same  as  that  described  for  NaCl  in  the  previous  chapter. 
No  ci'ystals  were  grown  for  this  part  of  the  experiment,  however, 
but  instead  large  single  crystals  of  KCl  were  obtained  from  the 
Harshaw  Cheriical  Company,  which  v/ere  then  annealed  and  cleaved  much 
like  the  NaCl  crystals.  Neutron  irradiation  v/as  again  performed  at 
the  Argonne  National  Laboratory,  However,  since  the  relative 
abundance  of  the  potassium  isotope  which  was  activated,  namely 
is  less  than  7 (as  compared  to  1-00  Yq  for  Na^^  } , the  activities 
of  the  activated  KCl  were  considerably  smaller  than  those  for  NaCl, 
Upon  arrival  at  Urbana,  the  avei’age  activity  of  the  unshielded  KCl 
was  about  10  rar/hr  at  4 feet.  This  fact  and  also  the  shorter  half* 
life  of  K^-^  (12.44  hours)  placed  an  upper  limit  of  about  3 days  on 
the  diffusion  time.  As  a consequence,  no  diffusion  measurements 
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could  be  made  below  500*^C„ 

The  method  for  sealing  off  the  diffusion  specimens  was 
identical  to  the  modified  procedure  described  in  Chapter  I {pp*  13 
and  14)  and  illustrated  in  F:igura  2,  In  contrast  to  NaCl,  however^ 
it  was  found  that  at  high  temperatxires  the  KCl  crystals  were 
exceedingly  sensitive  to  connamiiationo  This  effect  was  in  the 
form  of  a pink  layer  extending  some  distance  into  the  crystal  from 
all  sideSo  It  v/as  especially  pronounced  in  preliminary  experiments 
where  stainless  steel  capsules  were  used,  but  could  also  be  observed 
when  the  specimens  were  mere  .y  sealed  in  vycor  which  had  not  been 

■At 

coated  with  KGl  on  the  insid  The  KGl  coating,  the  hel^,  and  the 
platinum  capsules  inhibit  tho  reaction  between  the  specimens  used 
and  the  vycor-  w-it.h  all  the.se  precaution'’,  no  coloring  of  the 
crystals  could  be  observed  eren  when  left  at  high  temperatures  for 
extended  periods  of  timso  On  the  basis  of  these  facts  it  is 
believed  that  a metal  diffus  js  into  the  KGl  crystals  at  high  temper  - 
atureso  It  has  also  been  noted  that  the  pink  crystals  exhibited 
considerably  higher  self  dif  fusion  coefficients  than  the  clear  ones:, 
This  is  expected  if  di-  or  t;.-i-valent  metallic  impurity  ions  have 
diffused  into  the  lattice,  lending  further  strength  to  the  belief 
that  the  pink  color  is  due  t 3 such  ionSo 

5till  another  consequence  of  the  KCl  coating  of  the  vycor 
is  the  fact  that  the  samples  retained  the  fresh,  clear  appearance 
of  newly  cleaved  surfaceso  Samples  diffused  without  this  coating 
had  grayish,  dingy  looking  s u’faces,  and  in  particular  the  evapo- 
rated layftx"-  of  KGl  became  greyish  yellow  in  colors  Possibly  this 
is  due  to  .some  organic  matte':’  adhering  to  r.he  vycor 0 The  liquid 
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nitrogen  trap  also  helps  in  .his  respect,  perhaps  freezing  out 
organic  vapors  present  in  thtj  systemo 

The  slicing  of  the  KOI  specimens  was  carried  out  in  the 
same  manner  as  NaClo  However,  it  was  found  that  the  softer  KGl  was 
not  removed  as  uniformly  by  ';he  microtome  knife  as  NaCl,  leading  to 

O 

considerab?Le  scattering  of  data  on  the  log  activity  V3»  (depth) 
curves*  If  the  knife  orientation  across  the  crystal  surface  was 
along  the  110  direction  the  KGl  came  off  in  long  needles  rather  than 
in  powder  form,  and  effect iv?  slicing  was  almost  impossible*  But 
even  with  the  knife  oriented  in  a different  direction  the  slices 
had  to  be  taken  50  percent  t;iicker  than  for  NaCl*  Otherwise  there 
was  danger  of  the  scotch  tap  ; contacting  the  knife* 

The  errors  in  the  diffusion  coefficients  are  somewhat  large;* 
than  for  the  case  of  NaCl,  tlie  reason  being  primarily  the  increased 
thickness  of  the  slices.  This  has  the  effect  that  fewer  slices  can 
be  taken  (hence  larger  error  in  slope)  and  greater  approximations 
have  to  be  made  when  solving  the  diffusion  equation*  It  is  esti- 
mated that  the  errors  range  .'rcm  20  % at  700°G  to  40  % at  500°C* 
However,  these  errors  do  not  appear  sufficient  to  account  for  the 
scatter  of  the  data  shown  in  Figure  18*  It  is  not  known  at  the 
present  time  why  KGl  data  ar a inferior  in  precision  to  NaGl  data 
which  are  taken  by  the  same  personnel  using  the  same  equipmento 
The  conductivity  also  was  measu’^ed  by  the  same  methods 
discussed  in  Chapter  I*  There  was  no  change  in  conductivity  in  the 
intrinsic  region  during  cn>  jycling  process.  In  the  structure 
sensitive  region,  however,  tie  conductivity  increased  eonsiderably 
during  the  second  half  of  a cycle.  This  effect  is  attributed  to 


contamination  of  the  sample.  It  was  much  reduced  when  the  Improved 
sample  holder  shown  in  Figure  3 was  put  into  operation.  The  pre- 
cision of  the  conductivity  measurements  is  the  same  as  for  NaCl<, 
However,  since  there  was  no  volume  effect  and  no  concentration 
gradient,  successive  measure/aents  on  different  samples  in  the 
intrinsic  region  were  within  a few  percent  of  one  another* 

Co  Results  and  Discussion: 

The  conductivity  measurements  made  on  foxir  different  speci- 
mens are  shown  in  Figire  17«  The  specimens  are  not  distiiguished 
in  this  graph  because  of  the  uniformity  of  behavior,  slignt  devi- 
ations occurring  only  below  4C0®C.  The  ”knee**  of  the  curve  is 
located  near  450*^0.  This  gives  an  indication  of  the  purity  of  the 
crystals.  The  concentration  of  vacancies  due  to  impurity  ions  is 

of  the  order  of  P.  x 10^^/cc,  This  number  has  been  obtalnad  by 

(29) 

taking  the  data  of  Kelting  and  Uitt  (who  plotted  the  concentratiaa 
of  thermal  vacancies  vs.  temperature)  and  by  assuming  bhai  at  the 
knee“  the  concentration  of  thermal  vacancies  is  equal  bo  the 
concentration  of  impurity  vacancies.  Harshaw  WaCl  crystals  used 
in  pi’evious  vrark  have  been  found  to  contain  an  impurity  vacancy 
content  of  so  that  the  crystals  used  in  the  present  work 

are  appreciably  purer. 

Figure  IH  shows  the  results  of  ei^teen  diffusion  experi- 
ments covering  the  temperature  range  from  520®C  to  660®C,  The 
solid  curve  represents  the  diffusion  coefficient  calculatad  from 
the  conductivity  by  means  of  the  Einstein  relation.  Since  the 
contribution  of  the  chloride  ion  to  the  conductivity  is  appreciable 
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in  KCl  at  high  toinpe natures,  the  transference  number  data  of 
/ “SO  J 

Kerkhoff'^  were  used  to  obtain  the  contribution  due  to  the 
potass iujn  ion  alone.  . . ' 


In  view 

of  this  contribution  of  the  chloride  ion,  we  must  use  only  the 
conductivity  due  to  the  K"*”  ion  in  the  Einstein  relation.  The 
total  conductivity  ^ (shown  in  Figure  17)  is  related  to  by 


O’. 

% - 

CT“ 


(2/J 


vdiere  t^  is  the  transference  number  of  the  K”*"  ion.  For  temperatures 
in  excess  of  620®C  the  calculated  diffusion  coefficient  is  repre- 
sented by  a dashed  line  since  it  was  necessary  to  extz’apclate 
Kerkhoff's  data  above  this  temperature. 

From  the  slope  of  the  line  in  Figure  IB  we  ca;.  get  the  sum 
of  the  activation  energy,  U,  and  the  energy  of  formation,  , of  a 

single  positive  ion  vacancy.  Unfoirtunately . neither  U ncr  ^ are 

7 

known  separately  with  great  accuracy.  Using  the  value  of  1.20  eV 

for  W from  Kelting  and  Uitt^^^^  and  l.Bb  eV  for  U 4 from  the 

present  work,  we  find  U to  be  .6B  eV« 

•Cn  this  connection  it  rai^t  be  of  interest  to  mention 

(31) 

briefly  a calculation  by  Di<^ne5,  who  obtained  a theoretical 
value  for  the  activation  energy  for  the  jump  of  a pair  in  KCl. 

In  order  that  a neutral  pair  may  diffuse  it  is  necessary  that  ions 
of  both  charges  Jump  into  the  corresponding  vacancies.  Since  the 
negative  ion  is  nuch  less  mobile  than  tne  positive  one,  tae  rate* 
determining  step  must  be  the  migration  of  the  negative  ioa.  Now 

the  missing  posit,  ive  ion  lowers  the  barrier  opposing  the  movemant 


r 
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of  the  negative  ion  vacancy,  hence  there  is  a drop  in  activation 
energy  of  a pair  as  compared  to  that  of  a single  negative  ion 
vacancy.  Using  the  method  of  Bom  and  ilayer^^^^  as  modified  bv 
Mott  and  Littleton,  Dienes  arrived  at  an  activation  energy  for 
the  pair  of  c375  eV,  This  value  may  be  somewhat  in  error,  but 
it  seems  likely  that  the  activation  energy  does  not  exceed.  .5  eV, 

It  is  thus  considerably  below  that  for  a single  positive  ion  vacancy 


It  therefore  appears  that  pairs  may  have  rather  Icrgo 
mobilities.  However,  the  results  presented  in  Figure  IB  show  that 
the  Einstein  relation  is  satisfied  for  KCl  in  the  temperature  range 
covered  by  these  experiments.  Thus  even  if  the  mobility  of  the 
pairs  s.  ould  be  large,  their  concentration  is  probably  not  sufficient 
to  enable  them  to  contribute  appreciably  to  the  diffusion,. 

The  data  of  are  also  shovna  in  Figure  IB.  The 

present  data  do  not  agree  with  those  of  Witt.  V/e  believe  Witt*s 


results  to  be  affected  by  an  uiiknown  constant  error » 

In  conclusion  we  may  say  that  the  present  data  lend  no 


support  to  the  theory  that  vacancy  pairs  exist  in  KCl  in  the 
intrinsic  region.  V/hile  the  existence  of  such  pairs  is  not  ruled 
out,  eitlier  their  mobility  or  more  probably  their  concenti’ation 
Cor  both)  are  too  small  to  influence  the  self  di.^Tusxon  coefficient 
of  the  potassiinn  ion. 
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APPEi^DIX  I 

THE  EINSTEIN  RELATION 

We  shall  derive  the  Einstein  relation  first  by  a general 
thermodynamic  argument  and  then  from  an  atomistic  point  of  view 
vd.th  special  reference  to  diffusion  and  conduction  as  they  are 
supposed  to  occur  in  alkali  halide  crystals.  Since  the  experimen- 
tally measured  diffusion  and  conduction  are  essentially  one- 
dimensional, we  shall  limit  our  derivations  to  the  one-dimensional 
caser 

Let  us  suppose  we  have  particles  of  charge  e in  an  electric 
field  ^ ^ “ Then  the  density  of  particles  at  a point  x is 

given  by  Boltzmann’s  equation  as 

n(x)  = const  . (25) 

The  current  density  of  charged  pairticles  under  the  influence  of 
the  field  can  be  expressed  by 


1 = neuE  =“  -neu 


(26) 


where  u is  the  mobility.  This  current  flow  will  create  a concen- 
tration gradient,  and  thus  give  rise  to  a diffusion  current  opposing 
the  buildup  of  this  concentration  gradient.  The  diffusion  current 
density  is  given  by 


i - -e  D 


6n 

33E 


where  D is  the  diffusion  coefficient, 
current  is  flov/ing,  and  we  have 


(27) 


In  the  stationary  state  no 


-neu  4^  -e  D ^ » 0 
^ 3x 


(28) 
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By  integratii.g  (ZB)  we  obtain  immediately 


n = const  e a 
Comparing  (25)  and  (29),  we  arrive  at 


= e 

i?r 


(29) 


(30) 


which  is  the  form  the  Einstein  relation  is  often  given.  To  obtain 
the  form  used  in  the  present  work,  we  substitute  d*  “ neu  in  (30). 
We  thus  obtain 


cr  = ne 

IT 


(31) 


which  is  identical  to  (1). 

Next  let  us  consider  a more  detailed  model  consisting  of  an 
atomic  lattice  in  which  the  atoms  are  separated  from  their  nearest 
neighbors  by  a distance  a.  The  atoms  in  their  equilibrium  positions 
lie  in  energy  troughs  of  height  U.  The  probability  per  second  of  an 
atom  making  a jump  into  a neighboring  position  in  the  absence  of  a 
field  is  given  by 

V -U/kT  . 

P =■  4-  ~ e ' , (32) 

3 n 

In  this  expression  is  the  probability  that  the  atom  has 

enough  energy  to  overcome  the  barrier,  v is  the  vibration  frequency 

of  the  atom,  i is  the  probability  that  the  atoms  will  be  vibrating 

in  the  x direction,  and  ill  gives  the  probability  that  the  adjoining 

n 

lattice  site  will  be  vacant.  Here  n^  is  the  density  of  vacancies 
and  n the  density  of  a toms . 

Now  let  us  assume  that  a concentration  gradient  exists  (e.g. 
some  of  the  atoms  are  radioactive).  The  flux  P of  such  " tagged” 
3.ijOinC!  S,  CI'O  3 o d plane  A normal  to  the  x direction  and  midway  between 
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two  atomic  planes  is  given  by  the  flux  of  such  atoms  jumping  across 

it  from  one  direction  minus  the  flux  of  such  atoms  jumping  the 

opposite  way*  Each  flux  is  given  by  the  number  of  tagged  atoms 
2 

per  cm  on  the  atomic  plane  adjoining  the  fictitious  plane  A times 
the  jump  probability  per  second  across  the  plane#  Hence 


P - an,* 


an. 


(33) 


‘2  P2I 

\fdiere  n^^*  and  n2*  are  the  densities  of  tagged  atoms  on  either  side 
of  plane  A,  and  p^2  arid  respective  jump  probabilities 

per  second#  In  the  absence  of  a field,  p^j^r  » - p.  However,  due 

to  the  concentration  gradient  n^*  ^ ^2**  Substituting  (32)  in  (33) 
we  obtain 

^ ^ = _a^£TY*aH  n’t*  V « ' (34) 


P - a2 


-U/kT  2 * n„  i -U/kT 

JC  ''  e = -a^grad  n*  — i 2.  e 

n J n 3 


Using  the  definition  of  the  diffusion  coefficient 


we  get 


D 


P = -D  grad  n’' 


!lv  1 
n 3 


(35) 

(36) 


The  conductivity  may  be  derived  in  a similar  manner# 
Suppose  a uniform  field  E exists  in  the  x direction.  The  current 
density  across  plane  A can  now  be  written  as 


.k 

lit 


i = eP  » e an^  p^2  ~ «an2'*‘  Pg^ 


(37) 


We  are  no  longer  concerned  with  tagged  atoms,  since  both  tagged 
and  untagged  atoms  have  the  sarae  charge  e#  Thus  n^^*  = “ n# 

However,  the  field  E alters  the  height  of  the  barrier „ increasing 
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it  by  fa  in  one  direction  and  decreasing  it  by  the  same  amount  in 
the  opposite  direction.  Hence  p^^  ^ ^21"  J’^P  probabilities 

per  second  will  be 

n,  g _-(«/«■  - |||) 


'12 


n 3 ® 


i3^) 


n 


and 


■21 


^ n^  •/  -{U/kT  + f^) 

s — — e 


n 3 


(39) 


Substituting  (3S)  and  (39)  into  (37),  we  obtain 

-U/kT 


^ 2 slnh  ^ e 


i = e a iXy  ^ axrm  ^ 


(40) 


Under  normal  laboratory  conditions  eEa  « skT,  so  the  hyperbolic  sine 
may  be  replaced  by  its  argument.  We  then  have 


i = 


, e2a2n^  ^ -U/kT 


£ 


3kT 


Using  the  definition  for  conductivity 


i = o-  E 


we  thus  get 


<r  = 


2 2 
e a^ 


n„  J -U/kT 


e 


3kT 


Comparing  (36)  and  (43),  we  again  arrive  at  (31), 


(a) 


(42) 


(43) 


I 

* } 
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APPENDIX  II 

THE  DIFFUSION  EQUATION  AND  ITS  SOLUTION 
The  basic  law  for  one -dimensional  diffusion  of  particles  is 
known  as  Fick*s  first  law  and  can  be  stated  as 


P = - 


D 61 


(44) 


P is  the  flux  of  particles  (number  crossing  unit  area  perpendicular 

to  the  X direction  per  unit  time),  D is  the  diffusion  coefficient, 

and  I is  the  density  of  particleso  (IVe  choose  the  symbol  I because 

it  represents  the  radioactivity  which  is  our  measure  of  the  density 

of  the  particles  as  it  applies  to  this  experimento) 

Equation  (44)  holds  for  the  stationary  state,,  From  it  we 

can  derive  another  lav/  which  will  apply  to  the  n on-stationary  state « 

Let  us  consider  two  planes  of  unit  area  at  positions  x and  x<*dx  and 

perpendicular  to  the  concentration  gradient  , We  can  write  the 

accumulation  of  particles  between  the  planes  as  F-  P dx  or 
6P 

“(■|^)dXo  As  a result,  the  concentration  increases,  and  we  can  also 
writs  the  accumulation  between  the  two  planes  as  (,^)dXo  Equating 
the  tv;o,  we  get 

If  D is  independent  of  the  concentration,  (44}  and  (45)  yield 


D - 61 
6x^  6t 


(46) 


This  is  known  as  Pick’s  second  lawo  It  is  this  equation  for  which 
we  must  find  a solutiono 

Our  experimental  arrangement  is  that  of  a radioactive  layer 
of  thickness  h deposited  on  one  end  of  a semi-infinite  parallelo- 


t 

I 


! 


I 


pipsd  extending;  into  the  positive  x direction.  No  diffusion  can 


occur  across  the  face  x = 0,  Further,  at  t = 0 all  the  radioactive 


atoms  are  concentrated  uniformly  in  the  radioactive  layer,  and  none 


are  present  in  the  rest  of  the  crystal.  Accordingly,  v/e  can  va*ite 


our  taitial  and  boundary  conditions  as 


l(x,0)  = \ for  h 


(47) 


0 for  X > h 


(43) 


61  (0,t)  - 0 


Equation  (46)  can  be  solved  by  means  of  the  standard  method 


of  seoaratlon  cf  variables,  Ihc  solution  can  be  written  in  the  forr. 


I(x,t)  = l^c^(a)  cosa  x ->  €2(0.)  sin  a x 


-a'^Dt 


In  view  of  condition  (48)  it  can  be  seen  at  once  that  C2(a)«0,  The 


most  general  solution  will  be  a superposition  of  solutions  of  the 


form  (49)  over  all  positive  veiucs  of  a.  It  can  be  written  as 


I(x,t)  = J c^(a)cos  ctx  da 

0 


(50) 


We  next  calculate  c-,  (a),.  Stai'ting  with 


l(x,0)  = C2(a.l  cos  0.x  da 
0 


(51) 


we  obtain,  using  the  Fourier  Int'jgral  Theorem, 

0^  K 

c^(a)=  ^ y c(x“  ,0)cosa:c*dx*  J ^52) 

0 0 


the  last  identity  following  from  (47)o  Thus  we  can  revnnite  (50)  as 


h r 00 


0 L 0 


ax*coso,x  e 


-a^Dt 


da  1 dx* 


(53) 
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The  integral  in  the  brackets  can  be  evaluated  using  the  trigometrie 
ident ity 

2cosax'cosox  =»  cosa(x*-x)  -f  cosa(x'^+x)  (54) 

and  Formula  506  in  Peirce’s  Integral  Tables o As  a result, 

(53)  becomes 

I(x,t)  -ia j e dx'  (56) 

2>T«I5t 

x-h 

While  (56)  is  che  exact  solution  to  the  diffusion  equation 
(46),  it  is  not  readily  suitable  for  evaluation  of  D from  the  exper-^ 
imental  data„  It  will  be  recalled  that  such  data  are  in  the  form 
of  log  radioactivity  vso  ( depth )'^  curves o If  the  thickness  h of  the 
original  radioactive  layer  is  small  comparec.  to  the  total  depth  of 
diffusion,  we  can  make  the  approximation 


x+h 

/ 


4Dt 


x-h 


d>:» 


/V 


(57) 


The  error  jjitroduced  thi’ough  this  approximation  will  be  discussed 
below.  In  view  of  (56)  and  (57)  we  can  write 


Inl  ~ In  hl^  " _ In  hlo  (nX)^  (56) 

’(»Dt  40t  V«Dt'  4Dt 

where  X is  the  thickness  of  ore  slice  and  n the  depth  in  the  cjrystal 

in  units  of  X«  By  choosing  two  slices  a and  b in  the  crystal,  at 

deofchs  n„  and  n,  and  with  activities  I and  I.  , we  can  derive  an 
e.  t a D 

expression  for  D from  (56),  It  is  given  by 

n - i-  GtLr»ila!_  .2  (59) 

4b  Inlg-lnl^j 


•n 


4S 

which  is  identical  Wxth  (19) o This  expression  is  suitable  for  the 
calculation  of  D from  the  experiiaental  data. 

We  must  now  consider  more  closely  the  approximations  made 
by  using  eq  (57)«  Generally  speaking,  this  equation  means  approxi- 
mating an  area  under  a Gaussian  curve  by  a rectangle,  both  being 
centered  about  the  same  ordinate  and  both  having  width  h.  It  can 


be  seen  at  once  that  the  smaller  the  value  of  h,  the  better  the 
approximation  will  be.  For  h =»  0 the  exact  solution  equals  the 
approximate  solution.  However,  the  product  Dt  also  is  of  importance.. 
If  we  regard  the  Gaussian  cux*ve  as  invariant  the  interval  of  inte- 
Sration  will  contract  as  Dt  increases  and  vice  versa.  Thus  the 
larger  Dt  the  better  the  approximation  will  be. 

Another  source  of  error  in  D is  the  finite  thickness  of  the 
slices.  For  the  purpose  of  a log  activity  vs.  (depth plot  it  would 
be  reasonable  bo  consider  the  radioactivity  I of  a slice  to  be  con- 
centrated at  or  near  the  center  of  su  ch  a slice.  However,  in  our 
plot  of  log  I vs.  n , n represents  the  depth  measured  always  at  the 
end  of  a slice.  This  error  can  be  minimized  by  plotting  log  I vs. 
(n-l/2)^  instead  of  log  I vs.  n^„ 

The  approximat ions  and  errors  just  discussed  are  shown 
graphically  in  Figures  15  and  15.  Each  Figure  represents  a different 
value  of  Dt.  Curves  1 and  2 represent  the  logarithm  of  the  exact 
solution  for  two  different  values  of  h,  plotted  against  To 

obtain  numerical  values  for  these  solutions,  they  were  written  in 
the  form 


X4h 


o 


2'hi'ih- 


( 


4Dt 


dx 


r h+x  „ h~xl 


(60) 


» I 
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They  could  then  be  evaluated  with  the  help  of  error  function  tables. 
Curves  3 and  4 show  the  logaritlui  of  the  approximate  solution  (as 
given  by  the  right  hand  side  of  (57)  ) plotted  against  and  (n-l/2 ) 
respectivelyo 

The  units  used  for  I in  Figures  I5  and  I6  are  arbitrary. 

Ail  constants  common  bo  both  exact  and  approximate  solutions  have 
been  neglected.  Since  we  are  only  interested  in  the  slopes  of  the 
curves  and  not  the  absolute  magnitude  of  log;  I,  some  of  the  curves 
have  also  been  shifted  up  or  dovm  arbitrarily  so  as  to  allow  the 
easiest  comparison.  The  pertinent  expressions  and  numerical  values 
used  5jn  Figures  15  and  16  are  listed  in  Table  III. 


Tab3 

e III 

Figure  15 

Db 

h 

X 

X 

n 

Curve 

1 

6x10“'^  cra'^ 

10"'^  cm 

5xl0“^cm 

nX 

integer 

Curve 

2 

6xl0“^cm^ 

5xl0'“^cm 

5xl0“’^cm 

nX 

integer 

Curve 

3 

6x10"^'' cra^ 

0 

5 xlO^^cm 

n>. 

integer 

Cuiure 

4 

6x10“^ cm^ 

0 

5xl0"^cm 

(n 

- 1/2 )X 

integer 

Figure  16 

n*. 

h 

X 

n 

Curve 

1 

2x10“’°  cm2 

10'"’'’  cm 

5x10”4ce 

nX 

integer 

Curve 

2 

2x10"“^^' cm2 

5xl0“'^''cra 

5xl0~^cm 

nX 

integer 

Curve 

3 

2x10“^ cm^ 

0 

5xlO~4cm 

nX 

integer 

Ciirve 

4 

2x!L0"’^cm2 

0 

5xl0“^cm 

(n 

“ 1/2  )X 

integer 

Inspection  of 

curve  1 in 

Figures  15 

and 

16  reveals 

1 a con- 

s id  enable  deviation  from  a straight  line  at  low  values  of  n^.  This 
deviation  is  due  to  i,he  relatively  large  value  of  h.  In  curve  2. 
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Avhere  h is  half  this  value ^ the  deviation  has  largely  disappeared. 
Comparing  either  1 or  2 with  either  3 or  4 we  notice  that  the  slopes 
of  the  exact  solution  are  always  steeper  than  those  of  the  approxi- 
nate  solution.  Similarly,  4 always  has  a steeper  slope  than  3. 
However,  all  the  curves  have  a tendency  to  straighten  out  and  be- 
come more  parallel  as  n^  increases.  This  is  to  be  expected,  since 
the  erixjr  due  to  h is  a surface  effect  and  since  the  relative 
difference  between  (n-l/2)  and  n will  decrease  with  increasing 
The  mnnerical  values  shorn  in  Table  III  correspond  to 
typical  experimental  situations.  To  give  an  idea  of  how  the  various 
approximation  affect  our  experimental  results,  Table  IV  lists  the 
percentage  errors  iji  D which  occur  due  to  these  approximations.  In 
each  case  D has  been  computed  froia  the  slopes  of  two  of  the  curves 
in  Figiires  15  and  16  and  the  percentage  difference  calculated. 


Approximat ion 

3 instead  of  1 

3 instead  of  2 

4 instead  of  1 
4 instead  of  2 


Table  IV 
^ error  in  D 

(Dt  = 6x10  ’^  cm^)  (Dt  = 2x10"^  cm^) 


1^ 

2^ 

1% 


It  should  be  pointed  out  that  the  slopes  for  the  case 
Dt  - 2x3.0"°  cm  were  measured  at  somewhat  lovmr  values  of  n than 
for  the  case  Dt  ==  6x10“^  cm^.  This  is  in  lirie  w:*th  experimental 
conditions,  since  the  former  case  gener'^lly  v;ill  not  yield  as  many 
measurable  tilices  as  the  latter.  Had  the  slopes  been  compared  at 
the  same  values  of  n^,  there  vrouLi  have  been  very  little  difference 
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Figure  Two-dinensional  representation  of  various  lattice  defects 
At,  Perfect  lattice 

Bo  Frenkel  defects  (vaccmcies  and  interstitial  ions) 


Schottky  defects  (positive  and  negative  ion  vacancies} 

Do  Divalent  ion  in  lattice  (divalent  positive  io;.'  replacing 
s-.onovalent  posit ion  and  creating  positive  ion 


vacancy } 


E,^ 


Fo 


Complex  in  lattice  (divaleit  positive  ion.  and 
ion  vacrmcy  bound  together) 

.paired  vaca.nc.ies  (positive  and  negative  icr.  va 
bound  to  a- e tvas  r ) 
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gure 


A,  o 


A o 

Go 

Do 


Experirr.9nr>3.1  arrangement  for  diffusion 

Alkali  halide  specimens 

Platinum  capsule  covering  specimens 

Evaporated  layer  of  radioactive  material  (exaggerated./ 
Platinum  spacer  to  separate  specimens  ( exa.ggerated ) 
Vycer  capusle  (filled  v.'fth  helium) 


Fo  Thin  coat5.ng  of  same  material  as  specimen 
Go  Capsule  hojeier  (thin- -v?alled  stainless  steel  tube,  closad 
off  at  bottoiT,  and  connect ?d  to  brass  disk  it  top) 

Ho  Chromel-Alumal  thenrocoupie 

For  ail  actual  diffusion  exper ;lment  the  ho.'der  is  inserted 
onto  a vacuum  tig^t  quartz  tube  running  tlirougit  a Hoskins 


furnace^ 
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Figure  3-'  Holder  for  conductivity  neasurenents 
An  Solid  platinum  electrodes 
3o  Vycor  tube  (6  mm  OD ) 

Cc,  Vycor  tube  {22  mm  OD) 

Do  Steel  3prlag3  (the  attacl.ed  platinum  wii'ea  constitute 
the  g:oound  lead) 

Eo  l-latinum-f-lri ti-nuin  Rhodium  (10  9o  i thermoccupLe 
Fo  Platinum  w-li'e  {high  lead) 

Go  Stupakoff  seal 
Ho  Amphe»’ol  socket 

lo  Sample  (actr.al  size?  1 inm  x 2 mm  :,z  2 mi-if  the  1 miri 

dimension  being  the  one  separat ijn g the  two  electr-cdes'f 
For  actual  conductivity  measurements  the  holder  ia  inserted 
into  a vacuum  tigl->.t  quartz  tube  rurmnig  through  a Hoskins 
furnace  and  pui’e  helium  circulated  tlu'oug):'.  the  tube  and 
holder c 
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Fifotre  4e 
Ac 
Bo 
Co 
Do 


Schematic  diagram  of  the  bridge  circuit 
H ewl ett-Pa ckar d aud io-o sc ill ator 
General  Radio  716»G  bridge 
Sample  and  sample  holder 
General  Radio  736-A  wave  analyser 


K-  liensral  Radio  7^^“C  precision  condenser 


All  the  condensers  and  resistors  shovm  inside  the  dotted 
Square  are  contained  within  the  bridge.  All  the  other 
equipment  is  connected  to  the  bridge  externally. 


i 


Flp;ure  Typical  plot  of  tha  logarifchia  of  the  radj  oactivity  I 
of  a slice  (in  ccuntxS  per  minute)  versus  the  square  of  the  slice 
number  n„  The  upward  deviation  of  the  first  few  points  is  due  to 
the  finite  thickness  of  the  original  radioactive  layer„ 
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Fifoire  5o  Degree  of  i'.Jjsociation  p calculated  as  a fimction  of 
T/Tq  for  single  crystals  of  NaCl  + ratio  ”*  30--'^-^ 

T * absolute  temperature 

Tq*  energy  of  association  (in  units  of  Boltzmann’s 
constant)  ’=  3900^^ 

Solid  curve;  p calculated  on  the  basis  of  dniation  (13-) 
Daslicd  cui'are;  p calculated  on  the  basis  of  equation  (7) 


Fif;ure  7q  Degree  of  association  p calculated  as  a fuiiCtion  of 
T/Tq  for  single  crystals  of  NaCl  + CdCl^  (molar  ratio  « I4xl0‘ 
T * absolute  ten.perature 


Tq=»  energy  of  association  (in  units  of  Boltzmann 
constant)  * 3900°K 


Solid  curve:  p calculated  on  the  basis 
Dashed  curve:  p calculated  on  the  basis 


of  equation  (11) 
of  equation  (7) 


Fifoire  Bo  Degree  of  association  p calculated  as  a function  of 
for  single  crystals  of  NaCl  CdClj>  (molar  ratio  ” 4x10  5) 

T » absolute  temperature 

Tq=«  energy  of  association  (in  units  of  Boltsroa.nn ’s 
constant]  “ 3900®K 

•Solid  ciiTves  p calculated  on  the  basis  of  equation  (11 ) 
Dashed  curve”  p calculated  on  the  basis  of  equation  (7) 


[ 
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Flp;ure  9»  Degree  of  association  p calculated  as  a function  of 
T/Tq  for  single  crystals  of  NaCl  + CdCl2  (molar  ratio  ""  ♦5x10 
T * absolute  temperatures 

energy  of  association  (in  units  of  Boltsraann’s 
constant)  ” 39CO°K 

Solid  curve:  p calculated  on  the  basis  of  equation  (11) 
Dashed  curve;  p calculated  on  the  oasis  of  equation  (7) 


Nad  + Cd0l2  S3  a function  of  the  inverse  abt.uiute  temperature^ 
for  various  molar  ratios  o 


Ao  Molar  ratio 
Be  Molar  ratio 
Co  Molar  ratio 

Do  Molar  ratio 
Eo  Molar  ratio 


Cd 

US 


Gd 

T7a 


v\JL 


Cd 

TIS 


» 30x10“^ 

« 14x10"' 

" 4x10  ^ 

- .5:cl0"' 

*s  -.»>  { « pure 


NaCl) 


Each  specimen  was  put  through  the  same  temperature  cyclOo 
The  upper  ciirve  for  each  crystal  refers  to  conductivities 
measured  during  the  first  half  of  the  cycle  (temperature 
increasing)  and  the  lov;er  curve  to  conductivities  measured 
during  the  second  half  of  the  cycle  (temperature  decreasing)# 
For  ^ pure®  NaCi  {E}  the  two  curves  coincide^ 
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Fifoire  llo  Calculated  and  measured  self  diffusion  coefficients  D 
(em^/sec)  of  Na^  in  w^ingle  crystals  of  NaCl  + CdCl^  (molar  ratio 
Cd/Ma  **  30x10°’^)  as  a function  of  the  inverse  absolute  temperature. 
To  calculate  diffusion  coefficients  from  conductivities  the  Einstein 
relation  was  used® 

Solid  curve;  self  diffusion  coefficient  calculated  from 
the  measured  corjductivity 

Dotted  curve:  self  diffusion  codfficient  calculated  from 

the  conductivity  resulting  if  there  were  no  association 
(p  calculated  using  equation  (?)  ) 

Dashed  curve:  self  diffusion  coefficient  calculated  from 

the  conductivity  resulting  if  there  were  no  association 
(p  calculated  using  equation  (11)  ) 

Circles;  self  diffusion  coefficient  measured  by  radioactive 


tracer 
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Figure  12„  Calculated  and  measured  self  diffusion  coefficients  D 
(cm^/sec)  of  Na^  in  single  crj'-stals  of  NaCl  + GdCl^  (molar  ratio 
Gd/Na  =■  lAxiO^^)  as  a function  of  the  inverse  ahaolu.te  temperaturoo 
To  calculate  diffusion  coefficients  from  conductivities  the  Einstein 
relation  v/as  usedo 

Solid  curve:  self  diffusion  coefficient  calculated  from 
the  measured  conductivity 

Dotted  curve;  self  diffusion  coefficient  calculated  from 

the  conductivity  resulting  if  there  were  no  association 
(p  calculated  using  equation  (7)  ) 

Dashed  curve:  self  diffusion  coefficient  calculated  from 

the  conductivity  resulting  if  there  were  no  association 
(p  calculated  using  equation  (11)  ) 

Circles;  self  diffusion  coefficient  measured  by  radioactive 
tracer 


D (cm^/Sec) 


1 


Fip:ui"e  14c  Calculated  and  measured  self  diffusion  coefficients  D 
(cm^/sec)  of  Na*  in  ainr^le  crystenls  of  NaCl  CdCi^  (radar  ratio 
Cd/Na  " c5x10"5j  as  a function  of  the  inverse  absolute  temperaturco 
To  calculate  diffusion  coefficients  from  conductivities  the  Einstein 
relation  was  usedo 


Solid  curve;  self  diffusion  coefficient  calculated  from 
the  measured  conductivity 

dashed  curve:  self  diffusion  coefficient  calculated 
from  the  conductivity  resulting  if  there  were  no  asso- 
ciation (p  calculated  using  either  eriuation  (7)  or  (11),) 
GirclesS  self  diffusion  coefficient  measured,  by  radioactive 
tracer. 
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Figure  15 » Solutions  of  the  diffusion  equation 

Curve  1:  logarithm  of  exact  solution  (h“10  cm)  plotted 

Curve  logarithm  of  exact  solution  (h^axlO^^cm)  plotted 

against  n^ 

Curve  3 1 logarithm  of  approximate  solution  (h=0)  plotted 
2 

against  n 

Curve  4:  logaritlim  of  approximate  solution  (h“0)  plotted 

2 

against  (n-l/2) 

h “ thickness  of  original  radioactive  layer 
n ” slice  number 

The  units  used  for  the  solution  of  the  diffusion  equation 
are  arbitrary.  Some  of  the  curves  have  been  shifted  up 
or  down  to  permit  easier  comparison  of  their  slopes.  All 
solutions  were  evaluated  for  Dt  * 6x10  ^ cm^,  where  D is 
the  diffusion  coefficient  in  cm^/sec  and  t the  diffusion 
time  in  seconds. 


/ 


S2 


Solutions  of  the  diffusion  equation 
Curve  1:  logaritlirc  of  exact  solution  (h-lO'^cm)  plotted 
against 

Curve  2:  logarithm  of  exact  solution  (h-5xl0**‘^cm)  plotted 

against 

Curve  3:  logarithm  of  approximate  solution  {h«0)  plotted 

against  n^ 

Curve  4;  logarithm  of  approximate  aolut^nr. 
against  (ri«l/?)~ 

h =•  thickness  of  original  radioactive  layer 
n » slice  number 

The  units  used  for  the  solution  of  the  diffusion  equation 
are  arbitrary.-.  Some  of  the  curves  have  been  shifted  up 
or  do;vn  to  permit  easier  comparison  of  their  slopes.  All 
solutions  were  evaluated  for  Dt  *•  2x10“^  cm^,  where  D is 
the  diffusion  coefficient  it:  cm^sec  and  t the  diffusion 
time  in  seconds. 


’yrnT' 


wry,  AKBfTfiAJiy  units) 
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Fif^re  17o  Conductivity  (ohm  cm  ) in  **  pure*  Harshaw  KCl 
sin£ie  crystals  as  a function  of  the  inverse  absolute  temperaturec 
The  circles  represent  the  results  cf  actual  ii.easureriients  on  four 
difrerent  specimens..  The  ■ knee®  of  the  curve  occurs  at  about 
450®Co 
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Flfxire  1^6  Calculated  and  measured  self  diffusion  coefficients 
D(cm'^/sec)  of  K"*"  in  ” pure”  Harshaw  KCl  single  crystals  as  a 
function  of  the  inverse  absolute  temperature.  To  calculate  self 
diffusion  coefl icierius  from  comluct iviticc  the  Einstein,  relation 

was  usedi, 

Solid  line:  self  diffusion  coefficient  calculated 
from  the  conductivity  due  to  the  potassium  ion 
only  (ur^ng  the  transference  number  data  of  Kerkhoff ) 
Dashed  r lf  diffusion  coefficient  calculated 

from  the  conductivity  due  to  potassium  ion  only 
(using  tho  extrapolated  transference  number  data 
of  Kerkhoff) 

Open  circles:  self  diffi-.slon  coefficients  measured  by 
radioactive  tracer  (present  v^ork) 

Dark  circles:  self  diffusion  coefficients  measured  by 


radioactive  tracer  (Witt) 
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